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Summary. In this paper we study the existence of stationary solutions for stochastic partial differential 
equations. We establish a new connection between Lp(R d ; R 1 )<g)I/ 2 ,(R d ; R d ) valued solutions of backward doubly 
stochastic differential equations (BDSDEs) on infinite horizon and the stationary solutions of the SPDEs. 
Moreover, we prove the existence and uniqueness of the solutions of BDSDEs on both finite and infinite horizons, 
so obtain the solutions of initial value problems and the stationary solutions (independent of any initial value) 
of SPDEs. The connection of the weak solutions of SPDEs and BDSDEs has independent interests in the areas 
of both SPDEs and BSDEs. 
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1 Introduction 

Let u : [0, oo ) x U x fl — > U be a measurable random dynamical system on a measurable space (U, B) 
over a metric dynamical system (J?, P, (Ot)t>o)i then a stationary solution is a & measurable 
random var iable Y : Q, -> U such that (Arnold 1 ) 

u(t,Y(u>),u)=Y(6tw) for all t > a.s.. (1.1) 

This "one-force, one-solution" setting is a natural extension of equilibria or fixed points in deterministic 
systems to stochastic counterparts. The simplest nontrivial example is the Ornstein-Uhlenbeck process 
defined by the stochastic differential equation du(t) = —u(t)dt + dB t . It defines a random dynamical 
system u(t,uo) = u^e^ 1 + j*e~^ t ~ s " > dB a and its stationary point is given by Y(u) = f® oo e s dB s . 
Moreover, for any uq, u(t, uq, 9-toj) — * Y(co) as t — > oo, where 9 t is the shift operator of the Brownian 
path: (9 t B)(s) — B(t + s) — B(s) for any s € (— oo,+oo). A pathwise stationary solution describes 
the pathwise invariance of the stationary solution over time along the measurable and P-preserving 
transformation 6t: fi — > fi, and the pathwise limit of the solutions of random dynamical systems. 
Needless to say, it is one of the fundamental questions of basic importance ([T], [7], [2], [2T], [50])- 
For random dynamical systems generated by stochastic partial differential equations (SPDEs), such 
random fixed points consist of infinitely many random moving invariant surfaces on the configuration 
space due to the random external force pumped to the system constantly. They are more realistic models 
than many deterministic models as it demonstrates some complicated phenomena such as turbulence. 
Their existence and stability are of great interests in both mathematics and physics. However, in 
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contrast to the deterministic dynamical systems, also due to the fact that the external random force 
exists at all time, the existence of stationary solutions of stochastic dynamical systems generated e.g. 
by stochastic differential equations (SDEs) or SPDEs, is a difficult and subtle problem. We would 
like to point out that there have been extensive works on stability and invariant manifolds of random 
dynamical systems, and researchers usually assume there is an invariant set (or a single point: a 
stationary solution or a fixed point, often assumed to be 0), then prove invariant manifolds and stability 
results at a point of the invariant set (Arnold [1] and references therein, Ruelle [28], Duan, Lu and 
Schaumulfuss 10J , Li and Lu |19j . Mohammed, Zhang and Zhao [21] to name but a few). But the 
invariant manifolds theory gives neither the existence results of the invariant set and the stationary 
solution nor a way to find them. In particular, for the existence of stationary solutions for SPDEs, 
results are only known in very few cases ([7], [14], [21], [29], [30]). In [29], [30], the stationary strong 
solution of the stochastic Burgers' equations with periodic or random forcing (C 3 in the space variable) 
was established by Sinai using the Hopf-Cole transformation. In [3T], the stationary solution of the 
stochastic evolution equations was identified as a solution of the corresponding integral equation up to 
time +oo and the existence was established for certain SPDEs by Mohammed, Zhang and Zhao. But 
the existence of solutions of such a stochastic integral equations in general is far from clear. 

The main purpose of this paper is to find the pathwise stationary solution of the following SPDE 

dv(t, x) = [5£v(t, x) + f(x, v(t, x), a* (x )Dv{t,x))]dt 

+g(x, v(t, x), a*(x)Dv(t, x))dB t , (1.2) 

without assumption that there is an invariant set. Here B is a two-sided cylindrical Brownian motion 
on a separable Hubert space Uq; Jz? is the infinitesimal generator of a diffusion process X\' x (solution 
of Eq. (|2~TTj) ) given by 

1 n ffi n f) 

with (ay (x)) = aa*(x). Eq. (jl.2[) is very general, especially the nonlinear functions / and g can 
include V« and the second order differential operator Jzf is allowed to be degenerate, while in most 
literature, g is not allowed to depend on V« or g only depends on Vu linearly (Da Prato and Zabzcyk 
[8], Krylov [16], Pardoux [23] ). As an intermediate step, the result of existence and uniqueness of the 
weak solutions of (|1.2[) , obtained by solving the corresponding backward doubly stochastic differential 
equations (BDSDEs), appears also new. The existence and uniqueness of such equations when g is 
independent of Vu or linearly dependent of Vu were studied by Da Prato and Zabzcyk [8], Krylov 
[16 . But we don't claim here our results on the existence and uniqueness for the types of SPDEs 
studied in [8] and [16] have superseded their previous results. 

Note that from the pathwise stationary solution obtained in this paper, we can construct an invari- 
ant measure for the skew product of the metric dynamical system and the random dynamical system. 
In this connection, we mention that in recent years, substantial results on the existence and uniqueness 
of invariant measures for SPDEs and weak convergence of the law of the solutions as time tends to 
infinity have been proved for many important SPDEs ([5], [6], [8], [12], [13] to name but a few). The 
invariant measure describes the invariance of a certain solution in law when time changes, therefore it 
is a stationary measure of the Markov transition probability. It is well known that an invariant mea- 
sure gives a stationary solution when it is a random Dirac measure. Although an invariant measure 
of a random dynamical system on R 1 gives a stationary solution, in general, this is not true unless 
one considers an extended probability space. However, considering the extended probability space, one 
essentially regards the random dynamical system as noise as well, so the dynamics is different. See [20] 
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for some examples of SDEs on K 1 and a perfect cocycle on S 1 having an invariant measure, but not 
a stationary solution. In fact, the stationary solution we study in this paper gives the support of the 
corresponding invariant measure, so reveals more detailed information than an invariant measure. 

In this paper, BDSDEs will be used as our tool to study stationary solutions of SPDEs. We will 
prove that the solutions of the corresponding infinite horizon BDSDEs give the desired stationary 
solutions of the SPDEs (|1.2|) . Backward stochastic differential equations (BSDEs) have been studied 
extensively in the last 16 years since the pioneering work of Pardoux and Peng [24] . The connection 
between BSDEs and quasilinear parabolic partial differential equations (PDEs) was discovered by 
Pardoux and Peng in [25] and Peng in [27]. The study of the connection of weak solutions of PDEs 
and BSDEs began in Barles and Lesigne [4]. The BDSDEs and their connections with the SPDEs were 
studied by Pardoux and Peng in [26] for the strong solutions, and by Bally and Matoussi in [3] for the 
weak solutions. On the other hand, the infinite horizon BSDE was first studied by Peng in [57] and 
it was shown that the corresponding PDE is a Poisson equation (elliptic equation). This was studied 
systematically by Pardoux in [55] . Notice that the solutions of the Poisson equations can be regarded 
as the stationary solutions of the parabolic PDEs. Deepening this idea, it would not be unreasonable 
to conjecture that the solutions of infinite horizon BDSDEs (if exists) be the stationary solutions of the 
corresponding SPDEs. Of course, we cannot write them as solutions of Poisson equations or stochastic 
Poisson equations like in the deterministic cases. However, it is very natural to describe the stationary 
solutions of SPDEs by the solutions of infinite horizon BDSDEs. In this sense, BDSDEs (or BSDEs) 
can be regarded as more general SPDEs (or PDEs). 

As far as we know, the connection of the pathwise stationary solutions of the SPDEs and infinite 
horizon BDSDEs we study in this paper is new (section 2). We believe this new method can be used 
to many SPDEs such as those with quadratic or polynomial growth nonlinear terms. We don't intend 
to include all these results in the present paper, but only study Lipschitz continuous nonlinear term to 
initiate this intrinsic method to the study of this basic problem in dynamics of SPDEs. We would like to 
point out that our BDSDE method depends on neither the continuity of the random dynamical system 
(continuity means u(t, •, w) : U — > U is a.s. continuous) nor on the method of the random attractors. 
The continuity problem for the SPDE (|1.2p with the nonlinear noise considered in this paper still 
remains open mainly due to the failure of Kolmogorov's continuity theorem in infinite dimensional 
setting as pointed out by some researchers (e.g. [10], [21]). 

One of the necessary intermediate steps is to study the BDSDEs on finite horizon and establish 
their connections with the weak solutions of SPDEs (Sections 3 and 4). Our method to study the 
L^R^jR 1 ) <g> L 2 p (R d ;R d ) valued solutions of BDSDEs on finite horizon was inspired by Bally and 
Matoussi's approach on the existence and uniqueness of solutions of BDSDEs with finite dimensional 
Brownian motions ([3]). But our results are stronger and our conditions are weaker. We will solve the 
BDSDEs driven by the cylindrical Brownian motion and nonlinear terms satisfying Lipschtz conditions 
in the space L 2 p (R d ; R 1 )®L 2 p (R d ; R d ). We obtain a unique solution (Y 4 '', Z 1 /) G S 2 '°([t, T];L 2 p (R d ; R 1 ))® 
M 2 <°([t,T}; L 2 p (R d ;R d )). The result F*>' e S 2fi ([t, T];L 2 (R d ; R 1 ), which plays an important role in 
solving the nonlinear BDSDEs and proving the connection with the weak solutions of SPDEs (also 
BSDEs and PDEs), was not obtained in [3J. The generalized equivalence of norm principle (Section 2), 
which is a simple extension of the equivalence of norm principle obtained by Kunita f[17j). Barles and 
Lesigne ([!]), Bally and Matoussi ([3]) to random functions, also plays an important role in the proofs 
of our results. We believe our results for finite time BDSDEs are new even for BSDEs. 

In section 5, we will solve the BDSDEs on infinite horizon and in section 6, we study continuity of 
the solution in order to ensure that it gives the perfect stationary solutions of the SPDEs. 
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2 The stationarity of the solutions of infinite horizon BDSDEs and 
stationary solutions of SPDEs 

On a probability space (f2,^,P), let (B t )t>o an d (W t )t>o be two mutually independent Q- Wiener 
process valued on U and a standard Brownian motion valued on M. d respectively. Here U is a separable 
Hilbert space with countable base {e^}^; Q € L(U) is a symmetric nonnegative trace class operator 

oo 

such that Qei = Ajej and A; < oo. It is well known that B has the following expansion (0): for 

»=i 

each t 

oo 



#,■(*) = — = < S t , e 3 >t7, j = 1,2, • 
A, 



where 



are mutually independent real- valued Brownian motion on (Q, J?, P) and the series (|2.1|) is convergent 
in i 2 (X2, J 5 ", P). Let TV denote the class of P-null sets of We define 

jr t)T 4 ^ T ^ &w \jM, for < t < T; 

&t = ^sfoo ® ^ V/ ^ for t > 0. 

Here for any process (r? t ) t > , J^ t = a{rj r -r) s ;Q<s<r< t}, = ,^ t , ^% = Vt>o ^t,T- 

Definition 2.1 Let § be a Hilbert space with norm \\ ■ ||§ and Borel a-field S". For K £ R + , we denote 
by M 2 '~ K ([0, oo); S) the set of ® £P j measurable random processes {4>(s)} s >o with values on § 
satisfying 

(i) 4>{s) : Q — > S is ^ s measurable for s > 0; 
^ P[/ o oo e-^||0( S )||Id s ]<oo. 

>lZso we denote by S 2 '~ ([0, oo); S) £/ie set of ® & / 'J?" measurable random processes {ip(s)} s >o 
with values on § satisfying 

(i) ip(s) : j? — > S is JF S measurable for s > and ip(-,cj) is continuous P-a.s.; 

(ii) P[sup s > e-^||^(s)|jl] <oo. 

Similarly for < £ < T < oo, we define M 2 '°([i,T];§) and S 2 '°([i, T]; §) on finite time interval. 

Definition 2.2 Let § be a Hilbert space with norm \\ ■ \\§ and Borel a-field . We denote by 
M 2 '°([t, T]; §) the set of £$n,T] ® & I measurable random processes {</j(s)}*<s<t with values on S 
satisfying 

(i) <j>(s) : fl — > § is ,^ S ,T V ^Too measurable for t < s < T ; 
(ii) E\jf \\(j){s)\\lds] < oo. 



v4/so we denote by S ' ([i, T];S) </«e set of S8[ t ,T] ® ^ / measurable random processes {ip(s)}t 
with values on § satisfying 

(i) i/j(s) : J? — > S is JF Sj t V ^Too measurable for t < s < T and tp(-, oj) is continuous P-a.s.; 
(ii) E[sup t<s<T yj(s)\\i] < oo. 



<s<T 
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For a positive K, we consider the following infinite horizon BDSDE with the infinite dimensional 
Brownian motion B as noise and Y t taking values on a separable Hilbert space H, Z t taking values on 
££ d (H) (the space of all Hilbert-Schmidt operators from R d to H with the Hilbert-Schmidt norm): 

/oc />oo 
e- Kr f(r 7 Y r ,Z r )dr + J K&- Kr Y r dr 

/OC p DO 

eT Kr g[r,Y r ,Z r )$B r - J e~ Kr Z r dW r , t > 0. (2.2) 

Assume / : [0, 00) x [2 x H x d (H) — ► H , g : [0, 00) X Q X H X C\ d (H) — ► C 2 Uq (H) are ® & ® 
8$h <8> im measurable such that for any (i, Yj Z) € [0, 00) x H x £jL (i?), /(t, Y, Z), g(t, Y, Z) are 

J^t measurable, where Uq — Q^(U) C U is a separable Hilbert space with the norm < u,v >u a = < 
Q~2u,Q~zv >u and the complete orthonormal base {vAie,-}?^, C^ o {H) is the space of all Hilbert- 
Schmidt operators from Uq to if with the Hilbert-Schmidt norm. It is noted that the Q- Wiener process 
{Bt)t>o is a cylinderical Wiener process on Uq, and both £f Io (H) and £ 2 , d (H) are Hilbert spaces. 

Note that the integral w.r.t. B is a "backward Ito's integral" and the integral w.r.t. W is a standard 
forward Ito's integral. The forward integrals in Hilbert space with respect to Q- Wiener processes were 
defined in Da Prato and Zabczyk [8]. To see the backward one, let {/i(s)} s >o be a stochastic process 
with values on C\j o (H) such that h(s) is & s measurable for any s > and locally square integrable, 
i.e. for any < a < b < 00, f b \\h(s)\\ 2 r2 ,fr\ds < 00 a.s.. Since & s is a backward filtration with 

respect to B, so from the one-dimensional backward Ito's integral and relation with forward integral, 
for < T < T, we have 

T fT'-t 

A, < h(s)ej, f k > Sp 3 {s) = - v/A" < h(T' - s)e h f k > d(3 3 {s), j, k = 1, 2, • • • 



where /%(s) = $j(T' — s) — (3j(T'), j = 1, 2, • ■ •, and so B s = Bt>- s — Bt>- Here {/&} is the complete 
orthonormal basis in H . From approximation theorem of the stochastic integral in Hilbert space ([8]), 
we have 

T'-t 00 j-T'-t 

h(T' - s)dB s = V / v/A~ < h(T' - s)ej, h > dPj{s)f k . 
v_ T j,k=i jT '- T 



Similarly we also have 

cT 



E 

j,k=l 



h(s)d^B s = V / v/A~ < ft(«) ei , / fc > rffaWh. 



It turns out that 



T /.T'-i 

h{s)d)B s = - h{T' - s)dB s a.s.. (2.3) 

Jt'-t 

Later we will consider another Hilbert space C^ o {H) (p > 2), a subspace of C 2 Io (H), including all 
h 6 C\j q (H) which satisfy 



h ^kw~ E ^\(hejj k )\ p < 00. 
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Definition 2.3 Let H be a dense subset of H . If (Y, Z) E S 2 >' K f| M 2 '~ K ([0, oo); H) (g) M 2 ^ K 
([0, oo); £^ d (H)), and for any p E H , 

/>oo />oo 

(e~ Kt Y u p) = ( e- Kr f(r, Y r , Z r )dr, p) + ( K e - Kr Y r dr, ip) 



e 



- Kr g{r,Y r ,Z r )SB T ,y)-{\ e~ Kr Z r dW r , ip), t > P - a.s., (2.4) 



or equivalently 

f (y tl vj) = (F T , ¥>) + (jf /(r, F r , Z r )dr, - ( jf 5 (r, F r , Z r )dt £ r , - (jf Z r dW r , , * 

\lim T ^ oo (e-^r T ^)=0 a.s., 

iften we call (Y, Z) a solution of Eg. [2.2$ in H. 

Remark 2.4 (i) Applying Itd's formula in H (see JSj). we have the equivalent form of Eq. 



Y t =Y T + jf f(r, Y r , Z r )dr - jf g{r, Y r , Z r )SB r - jf Z r dW r 
liniT^oo Q~ kt Yt = a.s.; 



(2.6) 



(ii) one can easily verify that the above definition doesn't depend on the choice of Hq due to the conti- 
nuity of the inner product; 

(Hi) the uniqueness of Y in S 2, ~ K ([Q, oo); H) implies if (Y',Z') is another solution, then Y s = Y£ for 
all s > a.s.. The uniqueness of Z implies Z s = Z' s for a.e. s S [0, oo) a.s.. But we can modify 
the Z at the measure zero exceptional set of s such that Z s = Z' g for all s > a.s.. 

The first main purpose of this section is to study the stationary property of the solution of BDSDE 
(|2.2| on H if the solution exists and is unique. In order to show the main idea, we first assume that 
there exists a unique solution of Eq. (|2.2p . The study of the existence and uniqueness of Eq. [|2.2|) will 
be deferred to later sections (sections 3-5). 

We now construct the measurable metric dynamical system through defining a measurable and 
measure-preserving shift. Let 8t ■ fl — ► ft, t > 0, be a measurable mapping on (Q, & ', P), defined by 
§t o B, = B s+t - B t , § t o W s = W s+t - W t . Then for any s,t>0, 

(i) P ■ 6- 1 = P; 

(ii) #o = /, where I is the identity transformation on ft; 

(iii) 9 S o § t = 6 s+t . 

Also for an arbitrary & measurable <p : Q — ► H , set 

eo(t>(uo) = <p(9(u)). 

We give the following bounded and stationary conditions for /, g w.r.t. 9.: 

(A.l). There exist a constant M x > 0, and functions /(•) € M 2 <~ A '([0, oo); R+), g(-) E M 2 <~ A '([0, oo); R+) 
s.t. for any s > 0, Y E H and Z E C 2 d {H), 



\\f{s, Y, Z)\\ 2 H < f 2 (s) + MiiiriH + Milieu 2 



\\g{s,Y,Z)\\U {H) <g 2 (s) + M 1 \\Y\\j 1 + M 1 \\Z\\ 2 c2 (H y, 

(A.2). For any r,s > 0, Y E H and Z E £^ d (H), § r o f(s,Y,Z) = f(s + r,Y,Z), 6 r o g{s,Y,Z) = 
g(s + r,Y,Z). 
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We start from the following general result about the stationarity of the solution of infinite horizon 
BDSDE. 

Proposition 2.5 Assume Ea A2.2\) has a unique solution (Y,Z), then under Conditions (A.l) and 
(A. 2), (Yt, Zt)t>o is a "perfect" stationary solution, i.e. 

9 r oY t = Y t+r , 6 r o Z t = Z t + r for all r, t > a.s.. 

Proof. Let B s — Bt'~ s — Bt> for arbitrary T' > and — oo < s < T' . Then B s is a Brownian motion 
with B = 0. For any r > 0, applying 9 r on B s , we have 

6 r ° B s — 9 r o (B T :_ S — B T i) = B T /_ s+r — B T , +r 

= {Bt'-s+t — Bt') — (Bx'+r — B T ) = B s ^ r — B- r . 

So for < t < T < T' and a locally square integrable process {h(s)} s >o, by 

r T'-t 



§ r O [ h( S )d)B s = -6 r O / h(T' - s)dB s 

Jt Jt'-t 

= - / 9 r o h{T' - s)dB s _ r 
Jt'-t 



IT'-T 
r T'-t— I 



pi — z — r 

I § r o h{T' - s - r)dB s 

JT'-T-r 
pT+r 

/ § r o h(s - r)<$B s . 

Jt+r 



It+r 

As T" can be chosen arbitrarily, so we can get for arbitrary T > 0, < t < T, r > 0, 

f-T r T+r 



pi pi +r 

§ r o / h(s)d i B s = § r oh(s- r)$B, 

Jt Jt+r 



(2.7) 



It is easy to see that <?(•, Y., Z.) is locally square integrable from Condition (A.l), hence by Condition 
(A.2) and (22) 

rT pT+r 

9 r o g(s,Y s ,Z s )dtB s = g(s,6 r oY s _ r ,§ r oZ s _ r )d}B s . (2.8) 

Jt Jt+r 

We consider the equivalent form Eq. (|2.6p instead of Eq. (12.2[) . Applying the operator r on both 
sides of Eq. (|2.6|) and by (|2.8[) . we know that 9 r o Y t satisfies the following equation 

9 r oY t = 6 r oY T + J^ r f(s, 9 r o y s _ r , § r o Z s _ r )ds 

~ / t +r 9(S, Or O Y s ^ r X o Z s _ r )tfB s - § r o Z s _ r dW s (2.9) 

Jim T ^ oo e- K ( T+ ^(9 r oY T ) = a.s.. 



On the other hand, from Eq. (|2.6j) . it follows that 

f Y t+r = Y T+r + f*+ r f(s, Y Sl Z s )ds - £+ r g(s, Y s , Z S )SB S - £+ r Z s dW s . , 

\ limT^oo e- K( - T+r ^Y T+r = a.s.. 1 ' ; 

Let Y. = 8 r o Y.- r , Z. = 6 r o Z.- r . By the uniqueness of solution of Eq. (J2T6J) and Remark 12.41 (iii), it 
follows from comparing (|2.9|) with (|2.10p that for any r > 0, 
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6 r oY t = Y t+r = Y t+r , 9 r o Z t — Z t+r — Z t+r for all t > a.s.. 

Then by perfection procedure ([T], [5]), we can prove above identities are true for all t, r > a.s.. We 
proved the desired result. o 

An important application of the BDSDEs is to connect its solution with the solution of the cor- 
responding SPDEs. If some kind of relationship is established, we can transfer stationary solutions 
from the infinite horizon BDSDEs to SPDEs. In this way, we are in access to stationary solutions of 
the SPDEs due to the stationary property of solutions of infinite horizon BDSDEs. For this, a specific 
Hilbert space H = L^R^R 1 ) defined below is considered. The main aim of rest of this section is 
to construct the stationary solution of the SPDEs. Some proofs are given in this sections. But many 
detailed proofs are postponed to later sections. 

In the following we consider the case H = L^R^R 1 ) with the inner product (7/1,1*2) = 
J Rd Ui(x)u2(x)p~ 1 (x)dx, a p-weighted L 2 space. Here p(x) — (1 + \x\) q , q > 3, is a weight function. It 
is easy to see that p(x) : R d — ► R 1 is a continuous positive function satisfying J Rd \x\ p p~ 1 (x)dx < 00 
for any p£ (2, q— 1). Note that we can consider more general p which satisfies the above condition 
and conditions in [3] and all the results of this paper still hold. 

We can write down the solution spaces following Definition 12.11 M 2 ~ K ([0, 00); L 2 (R d ; R 1 )), 
M 2 -- K '([0, 00); L 2 p (R d ;R d )) and S 2 ~ K ([0, oo);L 2 (R d ; R 1 )). Similar to the definition for M 2 ^ K 
{[0,oo);L 2 p (R d ;R d )) 7 we can also define M P '~ K ([0, oo);L p (R d ; R d )). 

For k > 0, we denote by Cf b (R p , R 9 ) the set of C fc -functions whose partial derivatives of order less 
than or equal to k are bounded and by Hp(R d ; R 1 ) the p- weighted Sobolev space (See e.g. 0]). In order 
to connect BDSDEs with SPDEs, the form of BDSDEs should be a kind of FBDSDEs (forward and 
backward doubly SDEs). So we first give the following forward SDE. 

For s > t, let Xl' x be a diffusion process given by the solution of 

X? = x + J* b(X^)du + £ a(X^)dW u , (2.11) 

where b G C ; 2 b (R d ;R d ), a £ Cf b (R d ; R d x R d ), and for < s < t, we regulate X f s > x = x. 
For any r > 0, s > t, x € R d , apply 6 r on SDE pTTT]) . then 

§ r o Xl' x = x + f + b(9 r o X^du + f + a(§ r o X l u x _ r )dW u . 

Jt+r Jt+r 

So by the uniqueness of the solution and a perfection procedure (c.f. [I]), we have 

§ r o Xl' x = XlXr X , for all r, s, t, x a.s.. (2.12) 

Moreover, it is well-known that the solution defines a stochastic flow of diffeomorphism Xp' : R d — > R d 
and denote by X\'' the inverse flow (See e.g. Kunita [17]). Denote by J(Xl ,x ) the determinant of 
the Jacobi matrix of X^ x . For (p S i/^R^jR 1 ), we define a process ip t : Q x [0,T] x R d -> R 1 
by if t {s,x) = ipiXl-^JiXl^). It is proved in [3J that ip t (s, ■) € F^R^jR 1 ) and for u e J ff^*(R d ;R 1 ), 
J Rd u(x)(p(x)dx = J2o<\ a \<k u a {x)D a tp{x)dx < J2o<\ a \<k \J ' Jr* \u a (x)\ 2 p- 1 (x)dxj Rd \D a ip{x)\ 2 p{x)dx < 
00 and f Kd u(y)<pt(s, y)dy = J Rd u(X l s ' x ) 
■ ip(x)dx. 

The following lemma plays an important role in the analysis in this article. It is an extension of 
equivalence of norm principle given in [18j . [J], [5] to the cases when tp and if - are random. 
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Lemma 2.6 (generalized equivalence of norm principle) Let p be the weight function defined at the 
beginning of this section and X be a diffusion process defined above. If s E [t, T], ip : fi x R d — > R 1 is 
independent of and (pp^ 1 E <g)R d ), then there exist two constants c > and C > such that 

cE[[ \^(x)\p- 1 (x)dx}<E[[ \ l p(X t s ' x )\p- 1 (x)dx]<CE[[ \<p{x)\p-\x)dx]. 



Moreover if $ : Q x [t, T] x R d -> R 1 , W{s, ■) is independent of and Vp- 1 E L X {Q <g> [t, T] <g> R d ), 
then 

cE[[ [ \<P{s,x)\p- 1 {x)dxds}<E[[ [ \^(s,X t s ' x )\p- 1 (x)dxds] 
Jt Js, d Jt Jm d 

<CE[[ [ \^(s,x)\p- 1 (x)dxds}. 



Proof. Using the conditional expectation w.r.t. and noting that - — p-^iy) " 1S ^<s measurable 
and \tp{y)\p~ l {y) is independent of &t s , we have 

E[[ \<p{Xl>*)\p-\x)dx] 



E[ E[ \rtv)\p-Hv) P ~ H ^/\*' V) \Ks] ]*V 

p (y) 

E\ W{v )\ p -\y) ]E [P-^p^^]dy. 

By Lemma 5.1 in [3], c < E[^^^^-} < C for any y E R d , s E [t, T], the claim follows. o 

By LemmaEH it is easy to deduce that X tr E M P '~ K ([0, oo); L^(R d ; R d )) for K E R+. 
Now we consider the following BDSDE with infinite dimensional noise on infinite horizon 

/oo />oo 
e- Kr f(X t r > x ,Y?< x ,Z t r > x )dr + J K G - Kr Y r ^ x dr 

poo poo 

-J e- Kr g(X£*,Yf> x ,Z*> x )d\B r - J e~ Kr (Z l r - x , dW r ). (2.13) 

Here B r = Y^jLi \/\ifij( r ) e j ; {/^j( r )}j=i,2,- - are mutually independent one-dimensional Brownian mo- 
tions. Note that we will solve Eq. (j2J3|) for F r *>' E L 2 p (R d ; R 1 ) and Z 1 / E C 2 ^ d {L 2 p {R d ; R 1 )) = L 2 p (R d ; R d ). 
Set gj = gy/Xjej : R d x R 1 x R d — ► R 1 , then Eq. ([2~T3|) is equivalent to 

/OO pOO 
e- Kr f{X t r x X ,X > Z r X )dr + J Ke- Kr Y r ^ x dr 

°° pOC pOC 

"E / e- Kr gj (Xp x ,Y r t>x ,K' x )^^(r)- e- Kr (Z^ x ,dW r ). 

Referring to Definition 1231 and noting that (^(R^R 1 ) is dense in L 2 (R d ; R 1 ) under the norm (/ Rd | • 
\ 2 p~ 1 (x)dx)i , we can define the solution in L^R^R 1 ) as follows: 

Definition 2.7 A pair of processes (Y.^^Z^) E S 2 -- R f| M 2 '- K ([0, oo); L^(R d ; R 1 )) M 2 -~ K 
([0,oo);L 2 p (R d ;R d )) is called a solution of Eq. flUSj) if for an arbitrary ip E C^R^R 1 ), 
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pOO p pOQ p 

e- Ks Y*' x tp(x)dx= / / e- Kr f(X t r ' x ,Y r t ' x ,Z t r ' x Mx)dxdr+ / / Ke- Kr Y*> x ip(x)dxdr 

J s JR d Js JM. d 

°° poo P 

J2 / e- Kr g J (X t r >*,Y^,Zp'Mx)dx<fl r 3 j (r) 

3 t 

( e- Kr Z^ x <p(x)dx,dW r ) P-a.s.. (2.14) 



Note that in (|2.14[) we leave out the weight function p in the inner product due to the arbitrariness of 
<P- 



If Eq. (|2.13[) has a unique solution, then for an arbitrary T, Y T ,X satisfies 



Y*' x = Y*' x + / /(X*' s , Y*> x , Z^ x )dr - / g(X}*, Y*> x , Z^d) B r - / (Zp x , dW r ). 

J s J s J s 

(2.15) 

In section 4, we will deduce the following SPDE associated with BDSDE (|2.15|) 



i(t, x) = u(T, x) + J [Jifu(s, x) + /(a;, u(s, x), (<r*Vu)(s, a;))]ds 

Sf(x, u(s, ar), (cr*Vu)(s, x))d^B s . (2.16) 



Here Jz? is given by (|1.3[) . u(T, a;) = l^'' 1 . But we can normally study general u(T, x) unless we consider 
the stationary solution. 

Now following Definition ^. 21 we write down the solution spaces needed in our paper: M 2 '°([t, T]; L 2 p {R d \ M 1 )), 
M 2 >°([t,T];L 2 p (R d ;R d )) and S 2 '°([t,T]; L^R^R 1 )). 

Definition 2.8 A process u is called a weak solution (solution in L 2 (U. d ;'U 1 )) of Ea. [KTB\) if (u,a*Vu) G 
M 2 '°([0,T]; L 2 p {M. d -M})) (g) M 2 >°([0,T]; L 2 (R d ;R d )) and for an arbitrary V G C7*'°°([0,T] x M^M 1 ), 

/ u(s, x)d s \P(s, x)dxds + / u(t,x)\f r (t,x)dx — / u(T,x)'l r (T,x)dx 

1 

2 

T 



T 

(a*Vu)(s,x)(cr*'7<l r )(s,x)dxds - I I u(s,x)div((b - A)W)(s,x)dxds 



f(x, u(s, x), (o-*'Vu)(s, x)})P(s, x)dxds 

oo r T 



-V) / / 9j(x,u(s, x), (o-*Vu)(s,x))>P(s,x)dxd' ! < Pj(s) P-a.s.. (2.17) 

This definition can be easily understood if we note the following integration by parts formula: for 

(pi,<P2€ C 2 (R d ), 

JiC<pi(x)(p2{z)dx — — I (a*V(fi)(x)(a*V(p2)(x)dx + / ipi(x)div((b — A)ip 2 )(x)dx. 

The main purpose of this section is to find the stationary solution of SPDE (|1 .2|) via the solution 
of BDSDE lpT3]) . We consider the following conditions: 
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(A.l)'. Functions / : R d x R 1 x R d — > R 1 and g : R d x R 1 x R d — ► Cfj^R 1 ) are 3§ Rd <g> <g> ^ R<i 
measurable, and there exist constants M 2 , M2J, C, Cj, ay > with Y^jLi < 00, Yl'jLi Q < 00 
and J2jLi a i < I s - t for an y y i' y 2 e iz^R^R 1 ), X 1 ,X 2 ,Z 1 ,Z 2 € £p(R d ; R d ), measurable J7 : 
R d ^[0,l], 

/ /7(x)|/(X 1 ( S ),r 1 (a;),Z 1 (. T ))-/(X 2 ( a ;),r 2 ( a ;),Z 2 ( a; ))| 2 p- 1 ( a ;)dx 
Jm d 

< [ U(x)(M 2 \X 1 (x) - X 2 {x)\ 2 + C\Y x (x) - Y 2 (x)\ 2 + C\Z x (x) - Z 2 (x)\ 2 ) p-\x)dx, 

JR d 

f U(x)\g j (X 1 (x),Y 1 (x),Z 1 (x)) - g j (X 2 (x),Y 2 (x),Z 2 (x))\ 2 p- 1 (x)dx 
Jm d 

< f U{x){M 2j \X 1 {x) - X 2 (x)\ 2 + CJY^x) - Y 2 (x)\ 2 + a^x) - Z 2 (x)\ 2 ) p~ l {x)dx; 

JR d 

(A.2)'. Forpe (2,g- 1), 

/ \f(x,0, 0)\ p p~ 1 (x)dx < 00 and / \\g(x, 0, 0)||% nl .p~ 1 (x)dx < 00; 

jR d jR d ^ ' 

(A.3)'. b € Cf 6 (R d ;R x ), cr € C ; 3 b (R d x R d ;R 1 ). Furthermore, for p is given in (A.2)', if L is the global 

Lipschitz constant for b and a, L satisfies K — pL — p ( p ~ x ) L? > 0; 
(A. 4)'. There exists a constant p> with 2p - pK - pC - p(p ~ 1) Cj > s.t. for any Yi, Y 2 € 

L^R^jR 1 ), A,Z € measurable C7 : R d -► [0,1], 

f U(x)(Y 1 (x) -Y 2 (x))(f(X(x),Y 1 (x),Z(x)) - f{X{x),Y 2 {x),Z{x)))p- l {x)dx 

JR d 

<-p [ U(x)\Y 1 {x)-Y 2 (x)\ 2 p- 1 (x)dx. 

JR d 

Remark 2.9 We need monotone condition (A. 4)' in order to solve the infinite horizon BDSDEs. But 
it does not seem obvious to replace the Lipschitz condition for f in (A.l)' by a weaker condition on f 
such as f is continuous in y using the infinite horizon BSDE procedure (e.g. [221). The difficulty is due 
to the fact that we consider various conditions in the space L 2 (M. d ; R 1 ) here rather than pointwise ones, 
therefore we cannot solve the BDSDEs pointwise in x. However, our conjecture is that the Lipschitz 
condition can be relaxed if we strengthen the monotone condition in Lp(R d ;R x ) to a pointwise one. We 
will study this generality in future publications. Here due to the length of the paper, we only consider 
the Lipschitz continuous function f to initiate this intrinsic method to the study of this basic problem. 

We first acknowledge the two theorems below and give their proofs in section 6. 

Theorem 2.10 Under Conditions (A.1)'-(A.4)', Ea. [2~Wjl has a unique solution (Y*' x , Z\' x ) . More- 
over £[su Ps > J* R- tr* K '\Y*>*\ p p- 1 {x)dx] < 00. 

Theorem 2.11 Under Conditions (A.1)'-(A.4)', let u(t,-) = Y t '' , where (Y'',Z.'') is the solution 
of Eg. H2.1S]) . Then for t £ [0, T], u(t,-) is a weak solution for Eg. &2.16\) . Moreover, u(t,-) is a.s. 
continuous w.r.t. t in Lp(R d ;R x ). 

Then we prove the main theorem in this section. 

Theorem 2.12 Under Conditions (A.1)'-(A.4)', let u(t,-) = Y t '' , where (Y.'',Z. , ') is the solution 
of Eg. £2.13]) . Then u(t,-) has an indistinguishable version which is a "perfect" stationary solution of 
Eg. (KM - 
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Proof. For Y S L^(R d ; R 1 ), Z e L^(K d ;K d ), let 

/(T, y, Z) = f(X*S,Y, Z), g(T, Y, Z) = g(X* a '\Y, Z). 
Here we take T = (s,t) as a dual time variable (t is fixed). By Condition (A.l)', we have 

||/(T,y,z)||i ?(Rd;R1) 

l/^Y^Z^lV 1 ^ 



<C P / \f(X t s x ,0,0)\'p- L (x)dx + C p \Y(x)\ z p- L (x)dx + C p \Z{x)\ z p-\x)dx. 

JR d JR d JR d 

Here and in the following, C p is a generic constant. By Lemma 121)1 and Condition (A. 2)', 

E[ / e- x *|/(^*0,0)| Vetted*] < C p / / e- K4, |/(rB,0,0)|V 1 (*)*B«ia 



/O JR d 

<C P [ \f(x,0,0)\ p p" 1 (x)dx < oo. 

JR d 

We take /(T) = (J Ri \f(X t s ' x ,0,0)\ 2 p- 1 (x)dx)i, then f(T,Y,Z) satisfies Condition (A.l). Similarly 
we can also prove g(T,Y,Z) satisfies Condition (A.l). On the other hand, applying 6 r on /(T, Y, Z), 
by flglj]) , we have for any Y e ^(K^K 1 ) and z £ L 2 p (R d ; R d ), 

K o /(T, y z) = f(9 r o X*--, y, z) = f(xiX r /, y z). 

Verifying g(T, y, Z) in the same way, we know that f(T, Y, Z) and g(T, Y, Z) satisfy Condition (A. 2). 
Since by Theorem l2.101 Eq. (|2.13p has a unique solution (Yj-, Z7-), this (Yr, Z7-) is a stationary solution 
as a consequence of Proposition 12.51 That is to say for any t > 



§ r oY T = 9 r o y*'' = y s *+ r '' , (9 r o Z T = 6>r ° = ' for a11 r > °) S > * a - S " 

In particular, for any £ > 

r o y/'' = y/+T'' for all r > a.s.. (2.18) 

By Theorem 12. 11[ we know that u(t, •) = Y t '' is the weak solution for Eq. (|2.16p . so we get from (|2.18|) 
that for any t > 

8 r o u(t, •) = u(t + r, •) for all r > a.s.. 

Until now, we know "crude" stationary property for u(t, •). And by Theorem l2.11[ u(t, •) is continuous 
w.r.t. t, So we can get an indistinguishable version of u(t, •), still denoted by u(t, •), s.t. 

9 r o u(t, •) = u(t + r, ■) for all t, r > a.s.. 

So we proved the desired result. o 



By Definition 12.81 Conditions (A.l)' and (A. 2)', one can calculate that g(-,u(s, •), (<r*Vu)(s, •)) 
G £f Jo (L'j l (R d ; R 1 )) is locally square integrable in [0, T]. Now we consider Eq. (11.2p with cylindrical 
Brownian motion B on Uq. For arbitrary T > 0, let Y be the solution of Eq. ()2.13|) and u(t, •) = Y t : ' 
be the stationary solution of Eq. (12.16|) with B chosen as the time reversal of B from time T, i.e. 
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B s = Bt-s — Bt or $j(s) — (3j(T — s) — (3j(T) for s > 0. By (|2.3p and integral transformation in 
Eq. (12.16j) . we can see that v(t, x) = u(T — t, x) satisfies (|1.2j) or its equivalent form 

v(t, x) = v(t, Vq)(x) = Vq(x) + / [Jifv(s, x) + f(x, v(s, x), (<7*Vf)(s, x))]ds 

Jo 

oo „ t 

+ V / g j (x,v(s,x),(a*Vv)(s,x))d/3 j (s), t > 0. (2.19) 

Here vo(x) — v(0, x). 

In fact, we can prove a claim that v(t, -)(uj) = Y^Zt''i^) does not depend on the choice of T. For 
this, we only need to show that for any T' > T, Y^Zt' '(&) — Y^Z*'' {&') when < t < T, where 
<2>(s) = B T _ S — B T and lj'(s) — B T '- S — B T >- Let 8. and 8[ be the shifts of &(■) and £/(•) respectively. 
Since by (j2. 18[> . we have 

y£3>-(&') = d' T ,_ t Y»'-{u') = Y^0 T ,_ t u>'). 

So we just need to assert that 0T-t& — 8 T ,_ t u)' . Indeed we have for any s > 

(§ T -tu)(s) = u>(T -t + s)- w(T - t) 

= {B T -(T-t+s) - B T ) - {B T ~(T~t) ~ Bt) 
= B t _ s — B t . 

Note that the right hand side of the above formula does not depend on T, therefore #T-t£>(s) = 
T ,_ t Lu'{s) = B t _ s - Bt. 

On probability space (f2,^,P), we define 0t = (0t) < , t > 0. Actually B is a two-sided Brownian 
motion, so (St)" 1 — 9-t is well defined (see pQ). It is easy to see that 6 t is a shift w.r.t. B satisfying 

(i) P-(0 t )- x =P; 

(ii) 9 = I; 

(iii) d s o6 t = 6s+u 

(iv) 9 t oB s = B s+t -B t . 

Since v(t, -){o-i) = u(T — t, ■)(£)) = Y^Zf ''0*0 a.s., so 

9 r v(t, -)(w) = <L r u(T - t, .)(<S) = u(T - 1 - r, ■)(&) = v(t + r, -)(w), 

for all r > and T > t + r a.s.. In particular, let = vq(uj) = Y™' (a)). Then above formula implies 

CUD: 

t F(a;) = F(6» t w) = w) = v(t, v (w),w) = v(t, Y(cj),uj), for all t > Q a.s.. 

That is to say 00*0 — ^(^t^XO — Y-r-t' 0*0 ^ s a stationary solution of Eq. (|1.2[) w.r.t. 6. Therefore 
we proved the following theorem 

Theorem 2.13 Under Conditions (A.1)'-(A.4)', for arbitrary T and t e [0,T], let v(t,-) = Y^ZlZ 
where (K*'', Z.'') is the solution of Eg. H2.13\) with B s = Bt- s — Bt for all s > 0. Then v(t, ■) is a 
"perfect" stationary solution of Eg. iTO)) . 
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3 Finite horizon BDSDEs 

Before we study the BDSDEs on infinite horizon, we need to study the BDSDEs on finite horizon 
and establish the connection with SPDEs. For finite dimensional noise and under Lipschitz condition 
for a.e. x £ R , the problem was studied in Bally and Matoussi [3]. In this section, we consider the 
following BDSDE with infinite dimensional noise on finite horizon: 

T i-T 

g{r, Xl' x , Y*< x , Zp x )d^B T - / (Z^ x , dW r ), < s < T. (3.1) 

J S 

Here h : Q x R d — ► R 1 , / : [0, T] x W l x R 1 x R d — > R 1 , g : [0, T] x R d x R 1 x R d — > C 2 Ua (R 1 ). Set 
9i ~ 9\/^ e 3 : [°> T ] x R d x R 1 x R d — ► R 1 , then Eq.^ is equivalent to 

Yy x = h(X^ x ) + F f(r, XP*, Y^ x , Zl' x )dr 

J s 

-J2 gj{r,Xp x ,Y}> x ,Zi> x )d)Ur)- {Zt' x ,dW r ), < s < T. 

j=l J s -Is 

We assume 

(H.l). Function h is ® SS^-d measurable and E[J Rd \h{x)\ 2 p~ 1 {x)dx] < oo; 

(H.2) . Functions / and g are ^[o,t] ®^R d 8><^ri ®2%w>- measurable and there exist constants C, Cj , ay > 
with J27LiCj < oo and oy < \ s.t. for any t G [0,T], Yi,F 2 € £p(R d ; M 1 ), X,Z X ,Z 2 G 

L 2 p (R d -R d ) 

f \f(t, X(x), Yi(x), Zi{x)) - f(t, X{x),Y 2 {x), Z 2 (x))\ 2 p-\x)dx 

JR d 

< C [ (\Yi(x) - Y 2 (x)\ 2 + \Zi(x) - Z 2 {x)\ 2 )p-\x)dx, 

JM d 

I \ 9j (t,X (x) , Yi (x) , Zi (x) ) - 9j (t, X (x) , Y 2 (x) , Z 2 (x))\ 2 p- 1 (x)dx 

JR d 

< f {C j \Yi(x)-Y 2 (x)\ 2 + a j \Zi(x)-Z 2 (x)\ 2 ) p - 1 (x)dx; 

JR d 

(H.3). J T / R d \ f(s,x,0,0)\ 2 p- 1 (x)dxds < oo and J Q T / Rd \\g(s, x, 0, 0)\\^ ^p--(x)dxds < oo; 
(H.4). b G Cf )fe (R rf ;R d ), cr G C ; 3 b (M d ;M d x M d ). 

Needless to say, the conditions (H.1)-(H.4) for the existence and uniqueness of solution of Eq. (|3.1[) 
are weaker than what are needed for the case of infinite horizon. We would like to point out that for 
the finite horizon problem, our conditions are weaker than those in Bally and Matoussi [3J. In (H.l), 
we allow the terminal function h depending on J^t.T independent sigma field oo- One can easily 
verify that it doesn't affect the results in [3]. Moreover, here we only need Lipschitz condition in the 
space L^R^R 1 ) instead of the pathwise Lipschitz condition posed in [3]. 

Definition 3.1 A pair of processes {Y. L ' , Z*.'') G S 2 '°([0, T]\ L 2 (R d ; R 1 )) M 2 >°([Q, T];L 2 p (R d ; R d )) is 
called a solution of Eg. S3. 1\) if for any <p G C° (R d ; R 1 ), 
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Y*> x ip{x)dx= I h(X t j, x )ip{x)dx + I I f(r,X$' x ,Y? ,x ,Zt.'*)<p(x)dxdr 



■ , Js JR d 



3=1 ' 
T 

Z*' x ip(x)dx, dW r ) P-a.s.. (3.2) 

The main objective of this section is to prove 

Theorem 3.2 Under Conditions (H.1)-(H.4), Eg. US. 1}) has a unique solution. 

This theorem is an extension of Theorem 3.1 in [3]. The idea is to start from Bally and Matoussi's 
results for finite dimensional noise and then take limit to obtain the solution for the case of infinite 
dimensional noise. But Bally and Matoussi's results cannot apply immediately here as we have a weaker 
Lipschitz condition and some of the key claims in the proof of Theorem 3.1 ([3]) are not obvious under 
iver, the result Y* 

a sequence of BDSDEs 



their conditions. Moreover, the result Yr G 5 2,0 ([0, T];L 2 (R d - 7 R 1 )) was not obtained in [3]. We study 



Y't.X.n Ufvt; x \ i / £ ( ~* \rt,X \rt.X.7l r/t,X,n\ t 

s ' = h(X T ) + / f{r,X r ' ,Y r ' ' ,Z r < )dr 

pT pT 



>>> pi pi 

-J2 Silr.^.^.^'V^r)- (Z^ n ,dW r ). (3.3) 

A solution of (031) is a pair of processes (Y t,v ", Z t, ' ,n ) G S 2 '°([0, T]\ L 2 p {R d ; K 1 )) (g) M 2 <°([0, T]; 
L 2 (IR !i ; E rf )) satisfying the spatial integral form of Eq. (l3.3j) . i.e. ()3.2|) with a finite number of one 
dimensional backward stochastic integrals. 
First we do some preparations. 

Lemma 3.3 Under Conditions (H.1)-(H.4), if there exists (Y.(-), Z.(-)) G M 2 '°([t, T];L 2 (R d ; R 1 )) 
^}M 2 '°([t,T];L 2 (R d ;W 1 )) satisfying the spatial integral form of Eg. EO)) for t < s < T, then Y(-) G 
S 2 '°([t, T]; L 2 (R d ; R 1 )) and therefore (Y s {x), Z s {x)) is a solution of Eg. EO)) . 



Proof. Let's first see Y s (-) is continuous w.r.t. s in L 2 p {R d ;R}). Since (Y s (x) . Z s (x)) satisfies the form 
of Eq.([3T5]) for t < s < T, a.e. x £ R d , therefore, 

\Y s+As (x)~Y s (x)\ 2 p- 1 (x)dx 



<C P I [ + \f(r,X t r ' x ,Y r (x),Z r (x))\ 2 drp- 1 (x)dx 



" p ps+As 

+CpE / I / g j (r,X t r ' x ,Y r (x) ) Z r (x))d^f3 j (r)\ 2 p- 1 (x)dx 

j =1 JR d Js 



s+As 



+C P / | / (Z r {x),dW r )\ 2 p- 1 (x)dx. 



For the forward stochastic integral part, it is trivial to see that for < As < T~ s, \ J^ +As (Z r (x), dW r ) \ 2 < 

su Po<As<T-s I J s S+Ai (Zr{%), dW r )\ 2 a.s.. And we can deduce that J Rd sup <A s <T-s | f^ +As (Z r {x), dW r ) \ 2 p~ 1 (x)dx < 
oo a.s. by the B-D-G inequality and Z.(-) G M 2 '°([i, T];L 2 p (R d ; R d )). So by the dominated convergence 
theorem, lim As ^ 0+ J Rd | f° + s (Z r (x), 
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dW r )\ 2 p~ 1 (x)dx = 0. Similarly we can prove lim As ^ - J Rd \ J^ +As (Z r (x), dW r ) \ 2 p~ 1 (x)dx = 
for t < s < T. The backward stochastic integral part tends to as As — > can be de- 
duced similarly. So Y s (-) is continuous w.r.t. s in L^R^R 1 ). From Conditions (H.2)-(H.4) and 
(Y(-),Z.(-)) £ M 2 ^°([t,T};L 2 (R d ;R 1 ))^)M 2 ' ([t,T}- 1 L 2 } (R d ;R d )), it follows that for a.e. x e R d , 

E[J^ r \f(r,Xy- c ,Y r (x),Z r (x))\ 2 dr} < oo and £J =1 E[fi \ 9j (r, Xp', Y r (x), Z r {x))\ 2 dr] < oo. For a.e. 
x E M. d , referring to Lemma 1.4 in [26], we use the generalized Ito's formula (c.f. Elworthy, Truman and 
Zhao [15]) to ip M (Y r (x)), where ip M {x) = x 2 I { _ M < x<M} + 2M(x - M)I {x > M} - 2M{x + M)I {x< _ M} . 

Til on 



Then 

tpu{Y s {x)) + / I{-M<Y r (x)<M}\Z r (x)\ 2 dr 



= ^ M (h(X^ x ))+j il;' M (Y r (x))f(r,X^,Y r (x),Z r (x))dr 

n -t 

+ J2 I{-M<Y r (x)<M}\gA^ X r X , Y r(x),Z r (x))\ 2 dr (3.4) 

3=1 Js 

n ,T *T 

~E / ^ M O^(x))g j (r,Xp 3) ,Y r (x),Z r (x))^0 j (r)- / ^ M (Y r (x))Z r (x),dW r ). 

We can use the Fubini theorem to perfect (|3.4|) so that (13. 4|) is satisfied for a.e. x € R d , on a full 
measure set that is independent of x. Taking integration in M. d on both sides, applying the stochastic 
Fubini theorem ([8j), we have 

tp M {Y s {x))p^ 1 (x)dx + / / I { ^ M <Y r (x)<M}\Z r {x)\ 2 p^ 1 (x)dxdr 



< 



iP M (h(X t T X ))p- 1 (x)dx+ / / ip' M (Y r (x))f(r,X t r : x ,0,0)p- 1 (x)dxdr 

Js JR d 

+ f [ ^M(Yr(x))(f(r,X t ^,Y r (x),Z r (x))-f(r,X t ^,0,0))p- 1 (x)dxdr 
Js Jm d 

+C P V f f \ gj (r,Xp x ,Y r (x),Z r (x)) - g j (r,X t r ' x ,0,0)\ 2 p- 1 (x)dxdr 
~{ Js JR d 

+ C pJ2 f I \9j(r,X t r ' x ,0,0)\ 2 p'\x)dxdr- [ ([ ^ M (Y r (x))Z r (x)p- 1 (x)dx,dWr) 
~z2 f i ^M(Yr(x))g j (r,X t r ' x ,Y r (x),Z r (x))p- 1 (x)dxd^ j (r). 



3=1 



Noting that ^.((M^ 1 )) < IM^t*)! 2 and \^'m( y t{x))\ 2 < 4|y r (a;)| 2 , so by Lemma [2H the B-D-G 
inequality and Cauchy-Schwartz inequality, we have 



E[ sup / i/j M (Y s (x))p 1 (x)dx] 

t<s<T jR d 

< C p E[ [ \h{x)\ 2 p- l {x)dx] + C p E[ [ [ {\Y r (x)\ 2 + \Z r {x)\ 2 )p- 1 {x)dxdr] 
./»<* Jt J« d 



+C P E{J2 / (\g J (r,x,0,0)\ 2 + \f(r,x,0,0)\ 2 )p- 1 (x)dxdr} 

] = 1 Jt JR d 
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-C P E\ 



\ 



\1>' M (Y a (x))\*p-i(x)dx / ]T \9j(r,X t r ' x ,Y r (x),Z r (x))\ 2 p-^x)dxdr] 



+C P EU [ [ \iP' u (Y s (x))\ 2 p-i(x)dx [ \Z r (x)\ 2 p- 1 (x)dxdr] 

V Jt JR d JR d 

pT r 



< C p £:[ / \h(x)\ 2 p- 1 (x)dx] + C p E[ / / (|y r (a;)| 2 + \Z r (x)\ 2 )p- 1 (x)dxdr] 

jR d Jt JR d 

n „T „ 

+C P E[Y] / / (\g J (r,x,0,0)\ 2 + \f(r,x,0,0)\ 2 )p- 1 (x)dxdr] 

j=1 Jt JR d 

+\e\ sup f \^ M (Y s (x))\ 2 p- 1 (x)dx}. (3.5) 

t<s<T jR d 

Since (K(-), £.(■)) € Af 2 '°([i, T]; L 2 (R d ; R 1 )) (g) M 2 '°([t, T];L 2 (R d ; R d )), taking the limit as M -> oo 
and applying the monotone convergence theorem, we have £'[sup t<s<T J Rd \Y s (x)\ 2 p _1 (x)dx] < oo. So 
Y.(-) e S 2 '°([i,T]; R 1 )) follows. That is to say (Y s (x), Z s {x)) is a solution of Eq.flO]). o 



For the rest of our paper, we will leave out the similar localization argument as in the proof of 
Lemma 13.31 when applying Ito's formula to save the space of this paper. 



Proposition 3.4 Under Conditions (H.1)-(H.4), assume Eg. 13.3\) have a unique solution (Yp x ' n , Zp x ' n ), 
then for any t < s < T , Yr' X °'*' n = Y*' x ' n and Zr' X '*' n = Z^ x ^ for any r E [s,T] and a.e. x E R d 
a.s.. 

Proof. For t < s < r < T, note that (Y" r s '''™, Zp' ,n ) is ^^oc ® measurable, so is independent of 
&Y S . Thus by Lemma we have 

E[[ T [ (\Y r s ' x '"'y + \Z s r ' x ^' n f)p- 1 (x)dxdr} 

J s JR d 

< C P E[ f [ (\Y r s ' x ' n \ 2 + \Z s r ' x ' n \ 2 )p~ 1 {x)dxdr] < oo. 

J s JR d 

Moreover, it is easy to see that Xp 3 = Xp x and (Y~r' s ' n ,Zp s '") is J^*, ® ^tr measurable, 

rS,Xt>',n „s,X* r ,n\ ^ „ /f 2.0^r„ TH. r2An>d. O. n^2.0/r„ tH. r 2 frnd. radw „„j fv s > x l »™ ryS^l'" ,n, 



so {Y. s ^°' ' n ,Z s ' A " e M 2 ' o ([s,r];i 2 ( R <i ;R i))0M 2 ^([ S ,T]; J L 2 (R <i ;R d )) and (F r s ' s '™,Z r s ^ s ' 
satisfies the spatial integral form of Eq.j33| for s < r < T. Define Yr' X '*' n = Yf' x > n , Z?' X °" = 
when t < r < s. Then (y r S:Xs Zr' Xs satisfies the spatial integral form of Eq. (|3.3p for t < r <T 
and (Y. a ' xt ''' n ,Zf' X '''' n ) E M 2 ' Q ([t,T};L 2 p (R d ;R 1 ))^)M 2 ' a ([t,T};L 2 p {R d ;R d )). Therefore, by Lemma 
(Yr' Xs '", Zr' Xs '") is the solution of Eq. (|3.3p . By the uniqueness of the solution of Eq. (|3.3p . we 



have for any s E [t,T], (Yr'*'*' 71 , Zr' X "'"' n ) = (F r *' x '™, Z^ x ' n ) for any r E [s,T] and a.e. x E R d a.s. 



Theorem 3.5 Under Conditions (H.1)-(H.4), Eg. i3.3\) has a unique solution, i.e. there exists a 
unique (K*''' n , Zf''' n ) E S 2 '°([t, T]; L 2 (R d ; R 1 )) M 2 '°([t, T]; L 2 (R d ; R d )) such that for an arbitrary 
if E C°(R d ;R 1 ) 

[ Y*' x ' n p(x)dx = [ h(X% x )ip(x)dx + f [ f{r,X t r ' x ,Y^' x ' n ,Z t r ^ n )ip{x)dxdr 

jR d jR d Js JR d 
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lp x > n <p(x)dx,dW r ) P-a.s.. (3.6) 

Proo/. Uniqueness . Assume there exists another (I*'*'", Z*> x > n ) € 5 2 '°([f, T]; L 2 p (R d ; R 1 )) M 2 <° 
([t,T];L 2 (R d ;R d )) satisfying Define F s *' x '™ = F s *' x '" - F s *' x '" and f*> x '" = Z*' x ' n - £*> x ' n , i < 

s<T. From Conditions (H.2)-(H.4) and {Y t ''' n , Zp'' n ) £ S 2 '°([t, T];L 2 p (R d ; R 1 )) Af 2 ^ 
([t,T];L 2 (R d ;R d )), it follows that for a.e. x € R d , P[/ 4 T |/(r, X*'*, r r *> x > n , Z*< x ' n ) - f(r,Xp x ,Y*< x < n , 
Zp x > n )\ 2 dr] < oo and £™ =1 \ gj {r,Xp x , Y^ n , Zp x < n ) - g 3 (r, Xp x , Y*< x ' n , Zp x > n )\ 2 dr] < oo. For 

a.e. x € R d , similar as in (|3.4|) . we use generalized Ito's formula to e Kr ipM(Y^' x,n ) where K £ R 1 , 
then take integration in W' x fl on both sides and apply the stochastic Fubini theorem. Note that the 
stochastic integrals are martingales, so taking the limit as M — > oo, we have 

E[e Ks f \ Y ^ n \ 2 p-\x)dx] +(K-2C-jr,Cj- \)E\ f f e Kr \Y^ n \ 2 p-\x)dxdr] 

JR d - =1 z Js JR d 

+(J-E^[/ T / e Kr \Zp x ' n \ 2 p- 1 (x)dxdr]<0. (3.7) 

z j=l Js JR d 

All the terms on the left hand side of (|3.7p are positive when taking K sufficiently large, so it is easy 
to see that for each s £ [t,T], Y^ ,x = a.e. x £ R d a.s.. By a "standard" argument taking s in the 
rational number space and noting J Rd e Ks \Y*' x ' n \ 2 p~ x (x)dx is continuous w.r.t. s, we have Y*' x > n = 
for all s £ [t,T], a.e. x £ R d a.s.. Also by (017]), for a.e. s £ [t,T], Zp x > n = a.e. x £ R d , a.s.. We can 
modify the values of Z at the measure zero exceptional set of s such that Zp x,n = for all s £ [t,T], 
a.e. x £ R d a.s.. 

Existence . Step 1: We prove for the following equation: 

i-T n rT pT 

Y^ n = h(X^ x )+ f(r,Xp x )dr-J2 9^,^)^$^)- (Zp x > n ,dW r ), (3.8) 

J S j — 1 

if (H.l) and (H.4) are satisfied, and f{;Xp'), gj(;X ty ) £ M 2 -°([t, T];L 2 (R d ; R 1 )), then there exists 
a unique solution. For this, we can first use a similar method as in the proof of Theorem 2.1 in [3] to 
prove there exists (Y?''' n , Zp^ n ) £ M 2 <°([t, T];L 2 (R d ; R 1 )) A/ 2 '°([£, T];L 2 (R d ; R d )) such that for an 
arbitrary tp £ C%(R d ;M}) 



a.s. 



Y*' x > n tp(x)dx = / h(X%*)(p(x)dx + / / f{r,Xp x )ip(x)dxdr 

JR d Js JR d 

n pT p pT p 

-V/ / gj(r,Xp x )<p{x)dx<#Pj(r) - / (/ Zp*< n <p(x)dx, dW r ) P 

j=l J s J Rd ">» J Rd 

By LcmmaGH Y t, ' ,n £ S 2 '°([t, T]; L 2 (R d ; R 1 )). Then Step 1 follows. 

Step 2: Given Qrt,»,n,iV-i j ^*,x J n,JV-ij e S 2 >°([t, T]; L 2 (R d ; R 1 )) M 2 -°([t, T]; L 2 p (R d ; R d )), define 
(y;*^-™^, Zp x ' n ' N ) as follows: 



1 / v t,x\ | / #»/ v - *^ \/-t.x.n.N — 1 ryt.x.n.N — 1\ ? 
*H T ) + y J( r , X r >V ; Z r ) dr 

-J2 g j (r,Xp x X' x,n,N ~\Z t r ,x,n ' N ~ 1 )Jp j {r)- (Zp x ^ N ,dW r ). (3.9) 

■ -i Js Js 
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Let (Y?' x > n >°,Zj.> x > n >°) = (0,0). By Conditions (H.l), (H.3), (H.4) and Lemma EH we know h, 
fir, Xp x , 0, 0) and gj(r, X*- x , 0, 0) satisfy the conditions in Step 1, so Eq. (|3.8p has a unique solution 
■ \: ,,A .Z ! - " A - g M 2 >%[t,T];L 2 p (R d ;R 1 ))®M 2 >°([t,T];L 2 p (R d ;R d )) when f(r, Xp x ) = f(r,X^ x , 0,0) 
and g{r,Xl' x ) — g(r, X*' x , 0, 0). From Proposition 13.41 and the Fubini theorem, we have Y^ x ' nA = 
yr,xp*,n,i and z t,x,n,i = z r,xp*,n,i for ag f g ^ x], x E R d a.s.. Thus by Conditions (H.l)- 

(H.4) and Lemma EH we have h, f(r,X*> x ,Y?' x > n < 1 ,Zp x ' n < 1 ) = f(r, Xp x , Yr' X *' m,n ' 1 , Z^**'*'"' 1 ) and 
g 3 [r, Xp x ,Y*> x ' n > 1 ,Z*' x > n ' 1 ) = ffj (r, X*> x , Y^**'"'"'^ Z^'"'"' 1 ) satisfy the conditions in Step 1. Follow- 
ing the same procedure, we obtain (Y. t ''' n ' 2 , Z*'''"' 2 ) e M 2 '°([t, T];L 2 (R d ; R 1 )) (g) M 2 >°{[t, T];L 2 (R d ; R d )). 
In general, we see flU) is an iterated mapping from S 2 '°([t, T]; L 2 p {R d ; R 1 )) (g) M 2 <°([i, T]; L 2 p (R d ; R d j) 
to itself and obtain a sequence {(Y*' x,n ' 1 , Zt' x ' n,l )}i = o,i.2---- We will prove that (|3 - 9(1 is a contraction 
mapping. For this, define 

\rt,x,n,i \/-t,x,n,i \rt,x,n,i — l r^t^x^n,i n,t,x,n,i >yt x^n^i — 1 

Y s — 1 s ~ 1 s' y A s — ~ A s ) 

f (s, x) = /(s, X*'*, Y*' x ' n,i , Z^<™< 4 ) - /(«, X* s ' x , Y*^- 1 , zl'^- 1 ), 

g){s,x) = ^J^n^Z*^ i = 1, 2, • • •, t < s < T. 

Then, for a.e. x e R d , (Y*' x ' n ' N , Zl' x ' n ' N ) satisfies 



nT _ n pT p r l 

Y^ N = f N -\r,x)dr-J2 gf- l {r,x)d)fc{r)- 

J S j—1 



(Zl' x ' n ' N ,dW r 



Applying generalized Ito's formula to e Kr \Y^ x ' n ' N \ 2 for a.e. x G Mr, by the Young inequality and 
Condition (H.2), we can deduce that 

[ e Ks \Y*' x ' n ' N \ 2 p- 1 {x)dx + K f [ e Kr \Y*' x ' n ' N \ 2 p- 1 {x)dxdr 

JR d Js JR d 

e Kr \Z t r ' x ' n ' N \ 2 p- 1 {x)dxdr 

e Kr (2C\Y r t ^ n -' N \ 2 + \\Y*' x ' n ' N - l \ 2 + \\Z t r ' x ' n ' N - l \ 2 )p- 1 {x)dxdr 
i 2 2 

e^^Qiy^-^ + ^^iz^- 1 ! 2 )^- 1 ^)^ 

~E r / e Xl '2y r *'^"^gf" 1 (r, a ;) /9 - 1 ( a ;)^t/3 J (r) 

-/ (/ e^y/^^^Z*' 31 '"^^- 1 ^)^,^). (3.10) 

Js JWL d 

Then we have 

(K - 2C)E[ f [ e Kr \Y^ x ' n ' N \ 2 p- 1 (x)dxdr] + E[ f [ e Kr \Z^ x ^ N \ 2 p- l {x)dxdr] 

Js J~R d Js JRd 

< ( = / e Kr ({l + 2Y,Cm- x - n ' N -\ 2 + \Zl- x ^ N ^\ 2 )p-\x)dxdrY 

Letting X = 1 + 2C + 2 Cj , we have 
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rT 



E[ / e Kr ((l + 2Y^C j )\Y r t < x < n - N \ 2 + \Z t r < x ^ N \ 2 )p- 1 (x)dxdr} (3.11) 
Js Jm d ~[ 

< i\ + f>)£[ f i e Kr ((l +2Y,C j m X ' 7hN - 1 \ % + \Z t r x - rhN - 1 \ 2 )p- 1 {x)dxdr]. 



3=1 j=i 

Note that ^[/ 4 T / Rd e^((l + 2X)~ t ' | 2 + I ' l 2 )^ 1 ^)^] is equivalent to E\jf j Rd (| • | 2 + 

I ' \ 2 )p~~ 1 {x)dxdr]. From the contraction principle, the mapping (|3.9p has a pair of fixed point 

(Y t '-' n , Z*'-'") that is the limit of the Cauchy sequence {(Y^ n ' N , Z^ n ' N )}^ =1 in A/ 2 >°([t, T];L 2 (R d ; M 1 )) (g) M 2 ' a ([t, 1 

We then prove Y t, '' n is also the limit of K*'''" ,JV in S 2 '°([t, T]; L 2 (R d ; R 1 )) as N -> 00. For this, we 

only need to prove {^ t ''' n ' iV }^ =1 is a Cauchy sequence in S 2 '°([t, T];L 2 (R. d ; K 1 )). Similar as in ([53]) . 

by the B-D-G inequality and Cauchy- Schwartz inequality, from (I3.10|) . we have 

E[ sup / e Ks \Y 6 t ' x ' n ' N \ 2 p~ 1 (x)dx] (3.12) 

t<s<T jRd 

< M 3 E[ f [ e^QS?'*'*'*- 1 ! 2 + + |^t,*.n,2V|2 + |^.".-V|2) p -l( a .) da . dr ] j 

where M3 > is independent of n and AT. Without losing any generality, assume that M > N. We 
can deduce from p. lip and (|3.12p that 

(E[ sup / |Y?' x ' n ' M -l^'*' n " w |V 1 (a!)d!B])* 

t<s<T jRd 
M . 

i=N+l t<*<TjR* 

< V (M 3 £[/ / e Xr (|y r t ' a ' n,i - 1 | 2 + |Z^' n ' i -Y + |y r t ' a: ' ,l,i | 2 + |Z*' a: ' n '' i | 2 )p- 1 (cc)d ; rdr]) 2 

< V (2M 3 E[ / e Kr ((l + 2^C J 0|y^ w -T + l^r"™' J "T)^ 1 (^)^) 5 

^ E b+E a i)^( 2il W / e^((l + 2^C J )|y r t ^ 1 | 2 + |Z^ 1 | 2 )p- 1 ( a; )^r])^ 
— ► as M, N — ► 00. 

The lemma is proved. o 

Following a similar procedure as in the proof of Lemma l3.3i and using Ito's formula to e Kr \Y^' x,n \ 2 , 
by the B-D-G inequality, we have the following estimation for the solution of Eq. (|3.3p : 

Proposition 3.6 Under the conditions of Theorem VS.'A (Y. '' ,n , Z. '''") satifies 

su P E[ sup f \Y*' x ' n \ 2 p- 1 (x)dx] + sup£[ f f \Z t r ' x ' n \ 2 p- 1 (x)dxdr] < 00. 

n t<s<T JM d n Jt JR d 

Remark 3.7 For s G [0,t], Eg. S3.3\) is equivalent to the following BDSDE 

yx,n = yt.x.n + j ^ ^ y *,„ Z *,n} dr 

/ gAr,x,r r X ' n ,Z x > n )d^ 3 (r)- / (Z x ' n ,dW r ). (3.13) 

■ -1 J s J s 
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Note that Y t ' x,n satisfies Condition (H.l). By a similar method as in the proof of Theorem \3. 5\ and 
Proposition^ we can obtain a (Y? n , Z? n ) e S 2 >°([0, t\; L 2 p {R d ; R 1 )) <g) M 2 '°([0, t\; L 2 {R d ; R d )), is the 
unique solution of Eg. [3.13]) . Moreover, 

sup £[ sup f \Y s x ' n \ 2 p- 1 (x)dx]+supE[[ f \Zp n \ 2 p~ l [x)dxdr] < oo. 

n 0<s<tJRd n Jo JR d 

To unify the notation, we define (Y*' x ' n ,Zl> x < n ) = (Y x ' n ,Z x ' n ) when s e [0,t). Then (Y*''' n , Z t '' ,n ) € 
S 2 '°([0,T]]L 2 (R d ;R 1 ))^M 2 '°([Q,T]-,L 2 p (R d ]M d )). Furthermore, we have 

sup E[ sup f \Y*' x ' n \ 2 p- 1 (x)dx]+swpE[[ [ \Z t r ' x ' n \ 2 p- 1 (x)dxdr] < oo. (3.14) 

n 0<s<T JR d n JO Jm d 



Proof of Theorem \3.2\ The proof of the uniqueness is rather similar to the uniqueness proof in Theorem 



Existence . By Theorem l3.5l and Remark l3.7[ for each n, there exists a unique solution (Y. Z t, ,n ) € 
^([O, T];L 2 (R d ; R 1 )) <g) M 2 <°([0, T];L 2 (R d ; R d )) to Eq.^. We will prove {Y t ''' n , £*'''") is a Cauchy 
sequence in S 2 '°([Q, T];L 2 {R d ; R 1 )) M 2 '°([0, T];L 2 (M. d ; R d )). Without losing any generality, assume 
that m > n, and define 



t,x,m,n \rt,x,m \rt.x.n ryt,x,m.n ryt.x.ra ryt,x,n 

i — X s ~ r s > — A s ~ A s i 

t,x \rt.x.ra r7t,x,m\ ff„ vt.x \rt,x,n nrt,x,n\ 



fm,n/„ \ -yt,x \rt.x.m rTt x.m\ r/ yt.x \^t,x,n ryt.x.n\ 

J \ s > x ) — j(s, a s , r s ,/j s ) — j{s,Jl s ,x s ,Zj s ■ ), 
9j {s,x) = 9j{s,X s > ,Y S ' ' ,Z S > )-g j (s,X s ' ,Y S > ' ,Z S ' ), < s < T. 

Then for < s < T and a.e. x £ R d , 

dY t,x,m,n = _/™,«( S; x ) ds + g™' n {a, x)d^ j [s) 

+ Z?= n+ i9As,Xi' x ,Y s t ' x > m ,Zt>*'™)dip j (s) + (Zt> x ^ n ,dW s ) 
Y T = U a.s.. 

Applying Ito's formula to e Kr ^Y^' x ' m ' n \ for a.e. x € R d , we have 

j e Ks\Y s t > x > m > n fp- 1 (x)dx + {K-2C- jrCj-h f [ e Kr \Y^ x ' m ' n fp- 1 (x)dxdr 
JR d - =1 2 J s J R d 

+(o-E a i)/ / z Kr \Zl^\ 2 p~\x)dxdr 

Z ~[ Js Jw d 

<C P / / {\Y^ m \ 2 + \Z^)p-\x)dxdr 

3 =n+l Js J* 1 

+ [ T [ \g ] (r 1 Xr,QM 2 p-\x)dxdr)}-Y j f [ 2e Kr Y^ n g^\r,x)p- 1 (x)dxd^ J (r) 

Js JK d j=l^ s ^ Rd 

rn -t „ 

- V / / 2e Kr Y r t > x > m < n g 3 (r,X^\Y r t > x > m ,Z t r < x ' m )p- 1 (x)dxd^ j (r) 

-Kr\rt-x.m,n >yt,x,m.n ~— 1 



(/ 2e Kr Y;' xm ' n Z t r ' x ' m - n p- 1 (x)dx,dW r ). (3.15) 

JR d 

All the terms on the left hand side of (|3.15j) are positive when taking K sufficiently large. Take 
expectation on both sides of (13. 15|) . then by Lemma [2761 and (|3.14p . we have 



22 Q. Zhang and H.Z. Zhao 

E[f [ e Kr \Y^ x > m - n \ 2 p- l {x)dxdr]+E[( [ e Kr \Z t r ' x ' m ' n \ 2 p- 1 (x)dxdr] 

JO JR d JO Jwi d 

m -T r 

<C V J2 ^ C 3 + ^)(^npE[ / / (|F r ^-"| 2 + \Z?' n \*)p-\x)dxdr] 

x, 0, Q)\ 2 p~ l {x)dxdr} } — > 0, as n, m — ► oo. (3.16) 
Also by the B-D-G inequality, from (|3.15|) we have 

E[ sup / e Ks \Y s t ' x ' m ' n fp- 1 (x)dx] 

0<s<T JK d 

<C P E[[ [ e Kr (\Y^' x ' m ' n \ 2 + \Z t r ' x ' m ' n \ 2 )p- 1 (x)dxdr] 

Jo Js. d 

+C P J2 (C j +a j )( Sa pE[[ T [ {\Y^> n \ 2 + \Z^\ 2 ) p -\x)dxdr\) 



>><■ />i p 

+C P V / \g j (r,x,0,0)\ 2 p- 1 (x)dxdr. 

, , , JO JR d 



j=n+l ' 

So by 1(333]) . (|3~16l) and Condition (H.3), we have 



E[ sup / e A ' s |Y g t ' a: ' m ' n | 2 /5- 1 (a;)da;] — > 0, as n, 

0<s<TjRi 



771 > OO. 



Therefore (K 4 '''™, Z.*'''™) is a Cauchy sequence in 5 2,0 ([0, T];L 2 (R d ; R 1 )) M 2 '°([0, T];L 2 (R d ; R d )) 
with its limit denoted by (Y^ X ,Z^ X ). We will show that (K 4 '', Z.*'') is the solution of Eq.([3"3|). i.e. 
(K*'', Z.*'') satisfies {3J]) for an arbitrary y> e C^R^R 1 ). For this, we will prove that Eq.(3j| converges 
to Eq. (l3.2|l in L 2 (fi) term by term as n — ► oo. Here we only show the convergence of the third term, 



j—lJs JR d 
oo „T p 

-J2 / 9 3 {r,Xl' x X' X ,Zl^{x)dxd)^{r)\ 2 ] 

J = 1 Js JR d 
n ,t f. 

]~~lJs JR d 

oo r T r 

+2E[\ J2 / / g J (r,X t r >*X'",%nrtx)<k<ft0i(r)\*] 
J=n+1 Js Js* 

oo r T r 

< c P £(Q + aj)E[ / / (\Y^ n - r r ^| 2 + \z^ n - z^fK 1 ^*-] 

~J JR d 

+C P E[ \ C I (g^X^X'^Kn-gj^X^^^dxd^^r)] 2 } 

J=n+1 Js Jm d 

+C V E[\ V / T / 3,(7,^,0,0)^)^^(^12]. 

„-_„_i_i Js JK d 



j=n+l 

Note 
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00 fT 



E l\ E / / (gj^x^X'^^-g^x^^^dxd^^r)] 2 } 

]=n+l Js Jr* 

= E[ 11/ (g^X^X'^Zpn-g^X^AO^tpWdxi £ Vi®^)*^*-] 
Js Jn* J=n+l 

= £[/ El / (^^'^^.^-S^X^O.O))^)^ V ^X]e 3 (e 3 ,e t )\ 2 dr} 

Js i=l "^ c! j=n+l 

= ^[E / T |/ (ft(r,^,^,^)- Sj (a;' i 1 o,o))^( a; ) ( fe| 2 *] 

00 „T 

< Q> E + f [ + \Z t r x \ 2 )R-\x)dxdr] — 0. (3.17) 

j=n+l " 

Here we used (X}^Ln+i Aj' e :7 ® e j) 3 = I^JLn+i \/\j e j ® e j- This can be verified as follows: for an 
arbitrary w e [/, by definition of tensor operator, 



( E V^jSj ® E V^ e J® e j) M = E V^7 e i( e j' E V^i e i( e h u )) 

j— n+1 i— n+1 j— n+1 i= n+1 



E v / ^7 e i(v / ^7 e J' e j)( e ^' u ) 



j=n+l 



= ( E Vj ® e j> 

j=n+l 



Similarly we have 



< 



C p^[I E / / 9 ] (r,X t r ^,0,0Mx)dxd^ J (r)\ 2 } 

j=n+l Js jRd 

V \g j (r,x,0,0)\ 2 p- 1 (x)dxdr — ► 0. (3.18) 



j'=ra+l 

That is to say (Y* ,x , Zl ,x ) Q< <T satisfies Eq. (|3.2p . The proof of Theorem 13.21 is completed. 



4 Weak solutions of the corresponding SPDEs 

In section 3, we proved the existence and uniqueness of the weak solution of BDSDE (|3.1|) . We obtained 
the solution (Y*< x , Z l / X ) by taking the limit of (Y^ x ' n , Z^ x ' n ) of the solutions of Eq.(0 in the space 
S 2 '°([0, T]; Lp(R d ; M 1 )) Af 2 >°([0, T}- 7 L 2 p (R d ; R d )). We still start from Eq.Q in this section. A direct 
application of Proposition [3]4] and Fubini theorem immediately leads to 



Proposition 4.1 Under Conditions (H.1)-(H.4), if we define u n (t, x) — Y t ' x,n , v n (t,x) — Zl ,x,n , then 



i n (s,X^ x ) = Y*> x > n , v n (s,Xl' x ) = Z^ x ' n for a.e. s £ [t,T], x e R d a.s. 



t 

d 
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We first use the idea of Bally and Matoussi [3] to give the correspondence between the weak 
solutions of SPDEs and BDSDEs with finite dimensional noise. Consider the BDSDEs (13.81). Define 



the mollifier K m (x) — mcexpj ^ rnx \y_ 1 }, if < x < K m (x) = otherwise, where c is chosen such 

that J+™K m (x)dx = 1. Define h m (x) = J Rd h(y)K m {x- y)dy, f m (r,x) = J Rd f(r,y)K m (x- y)dy and 
<7™(r, x) = J M<i gj(r, y)K m {x—y)dy. It is easy to see from standard results in analysis that h m {-) — > h(-), 
f m (r,-) /(ry) and 5™(ry) ^(ry) in L^M^M 1 ) respectively. Denote by (^&' n ,%i n ) the 
solution of the following BDSDEs: 

„T n pT pT 

Y?# B = h m (X% x ) + / f m (r, X*?*)dr - £ / ~gf{r, X^S^r) - / (Z%?, dW r ). 

J s j=l 

Let u m (t,x) — Y t ^ n . Then following classical results of Pardoux and Peng [26], we have Z t 'Jj n = 

<r*Vu£(t,a;), and = <(s,X^) = Y^'"'", = a* Vfi»(«, X*>*) - Z^""'". Moreover 

tt^(i, x) satisfies the smootherized SPDE. In particular, for any smooth test function <]/ £ C^ ,oo ([0, T] x 
R d ; R 1 ), we still have 

u% l (s,x)d s &(s,x)dxds+ [ (t,x)W r (t,x)dx - [ h m (x)&(T,x)dx 

JR d JR d 

\ I [ ((J*^u^){s,x)(<T*\7<lr)(s,x)dxds - [ [ u^(s,x)V({b-A)&)(s,x)dxds 

2 Jt JR d Jt JR d 

T p n pT p 

/ f m (s,x)&(s,x)dxds- V / / g^(s,x)^(s,x)dxd 1 $ j (s) P - a.s.. (4.1) 

JR d J = 1 Jt jR d 

But by standard estimates 

E[f I {\Y*% n - Yy x ' n \ 2 + \Z*>% n - Z t / x ^\ 2 ) P ~ 1 {x)dxds] — y as m -> oo. 



And as mi, to 2 — ► oo 

Jt JR d 

= E[f I (\Y^ - Y^\ 2 + \Z%£ - Zt^)p-\x)dxds) — » 0. (4.2) 

Jt jR d 

We define 7i to be the set of random fields {w(s, x); s £ [0, T], x £ M. d } such that (w, a*\7w) £ 
M 2 '°([0,T]]Ll(R d -R 1 ))<S)M 2 ' ([0,T];Ll(R d -,R d )) with the norm (E[j^ J Rd (\w(s,x)\ 2 
+ |((T*V)w(s, x)\ 2 )p~ 1 (x)dxds)i . Following a standard argument as in the proof of the complete- 
ness of the Sobolev spaces, we can prove H is complete. Now by the generalized equivalence of 
norm principle and ([4.2)1 . we can see that u 7 ^ is a Cauchy sequence in H. So there exists u n £ H 
such that ft,(T*Vfi:) -> (u n ,a*Vu n ) in Af 2,0 ([0, T];L 2 (W l ; R 1 )) M 2,0 ([0, T];L 2 (R d ; R d )). More- 
over Yf> x > n = u n (s,X*' x ), Zl' x ' n = a*Vu n (s,Xl' x ) for a.e. s £ [t,T], x £ R d a.s.. Now it is 
easy to pass the limit as to — > oo in (14.1)1 to conclude that u n is a weak solution of the cor- 
responding SPDEs. For the nonlinear case, we can regard f(r,x) — f(r,x,u n (r,x),a*S7u n (r,xj), 
gj{r,x) — gj{r,x,u n {r,x),a*Vu n (r,x)), and /, (jj satisfy the conditions in the above argument. Us- 
ing a similar proof as in the proof of Theorem 3.1 in [3] together with Theorem 13.51 and Proposi- 
tion SU we have, under Conditions (H.1)-(H.4), v n (t,x) = {a*Wu n )(t,x). Moreover, (u",cr*Vu") £ 
M 2 '°([O,T]; J L 2 (M d ;IR 1 ))0M 2 ' o ([O,r];L 2 (R' i ;R' i )), u n (t,x) is the weak solution of the following 
SPDE: 
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u n (t,x) = h(x) + J [&u n (s,x) + f(s,x,u n {s,x),(cr*\7u n )(s,x))}ds 

n „T 

-V / g j (s,x,u n {s,x),((j*Wu n )(s,x))d^ j (s), < t < s < T. 

That is to say, for any & G C^°°([0, T] x R d ; R 1 ), we have 

u n (s,x)d s <P{s,x)dxds + [ u n {t, x)&{t, x)dx - [ h(x)V(T,x)dx 

JR d JR d 

■\ I I {<T*Vu n ){s,x)(<j*W){s,x)dxds - [ [ u n (s,x)\7({b- A)W){s,x)dxds 
2 Jt JR d Jt Jm d 

f(s, x, u n (s, x), (cr*Vw")(s, x))W(s, x)dxds 



-V / / g j (s,x,u n (s,x),(a*Vu n )(s,x))^(s,x)dxd^P j (s) P - a.s.. (4.3) 
■ =1 Jt J& d 

By intuition if we define u(t, x) = Y t ,x , it should be a "weak solution" of the Eq. (|2.16p with u(T, x) — 
h(x). We will prove this result. 

First we need some necessary preparations. 

Proposition 4.2 Under Conditions (H.1)-(H.4), let (Y^ X ,Z^ X ) be the solution of Ea. Un)) . If we 

define u(t,x) = Y t ' x , then a*Vu(t,x) exists for a.e. t S [0, T], x g R d a.s., and u(s,Xl ,x ) = Yp x , 
(cr*Vu)(s,X*> x ) = Zl' x for a.e. s £ [t,T], x£R d a.s.. 

Proof. First we prove u n is a Cauchy sequence in TL. For this, by Lemma 12.61 and Proposition 14. 11 as 
m, n — ► oo , we have 

E[[ f (\u m (s,x) - u n {s,x)\ 2 + \(a*Vu m )(s,x) - {a*Vu n )(s, x^p- 1 (x)dxds] 

Jo JR d 

< C P E[[ T [ (K"( S ,X S <*)- U "( S ,X^)| 2 + |^ 

JO JR d 

= C P E[ [ [ QY°' x ' m - Y^ x > n \ 2 + \Z° S ' XM - Zf x < n \ 2 )p- l {x)dxds] — ► 0. 

Jo JR d 

So there exists u £ H as the limit of u n such that Vu(s,x) exists for a.e. s £ [0, T], x £ R d a.s. 
and E[Jq f Rd (\u n (s,x) - u(s,x)\ 2 + \(a* Vu")(s, x) - {a*Vu){s,x)\ 2 )p' 1 {x)dxds\ — ► 0. We define 
u(t, x) = Yt' x , then similar to the proof as in Proposition ^. U by the uniqueness of solution of Eq. (|3.1[) . 
we have it^X*' 31 ) = Y s Lx for a.e. s £ [t, T], x £ R d a.s.. Since 

E[ / \u(s, x) — u(s, x)\ 2 p~ 1 (x)dxds] 

Jo JR d 

<2E[[ [ (\u(s 1 x)-u n {s,x)\ 2 + \u n (s 1 x)-u{s,x)\ 2 )p- 1 (x)dxds] 

Jo JR d 

<C p E[f f (\Y S °' X -Y°> x > n \ 2 + \u n (s,x)-H{s,x)\ 2 )p- 1 {x)dxds} — -+ 0, 



u{t,x) = u{t,x) for a.e. t £ [0, T], x £ R d , a.s.. So a*\7u(t,x) exists for a.e. t £ [Q,T], x £ R d , a.s. 
Using Lemma 12.61 again, we have 
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rT 



E[f f (\u(s,X^)-Y^\ 2 + \(a*Vu)( S ,X^)-Z^\ 2 )p- 1 (x)dxds] 

Jt JR d 

<2E[[ [ Quis^l'*) - u n (s,Xl< x )\ 2 + \u n (s,X t s < x )-Y*> x \ 2 )p- 1 (x)dxds] 

Jt JR d 

+2E[ f T j (\(a*Vu)(s,Xr) ~ (a*Vu n )( S ,Xln\ 2 + \(<J*Vu n )(s, X^) - Z t - x \ 2 )p- l {x)dxds] 
Jt Jm d 

< C P E[ [ [ (\u(s,x) - u(s,x)\ 2 + \u{s,x) ~ u n (s,x)\ 2 + \Y^ x ' n - Y t ^ x \ 2 )p- 1 {x)dxds\ 
Jt Jm d 

+C p E[ [ [ (\(a*Vu)(s,x) - (a*Vu){s,x)\ 2 + |(<r*Vu)(«,x) - {a* Vu n )(a, x)\ 2 
Jt Jm d 

+\Z t s < x ' n - Z t s x \ 2 )p- 1 {x)dxds] — ► 0. 
So u{s,Xl' x ) = Y*' x , 0*Vu)(s,X*> x ) = Z* s ' x for a.e. s G [t,T], x G M d a.s.. o 

From Proposition ^. 21 and Lemma fZM it is easy to know that 
E[f [ \u n (s,x)~u{s,x)\ 2 p- 1 {x)dxds} 

Jt JR d 

+E[f [ \{a*\7u n )(s,x)-(cr*\7u)(s,x)\ 2 p- 1 (x)dxds} 
Jt Jm d 

<C P E[[ [ \u n (s,X t s ' x )-u(s,X t s ' x )\ 2 p- 1 (x)dxds] 
Jt Jm. d 

+C P E[[ [ \(a*Vu n ){s,X t s ' x )-(a*Vu){s,X t s ' x )\ 2 p- 1 (x)dxds] 
Jt Jwi d 

= C P E[[ [ \Y*' x < n -Y s t < x \ 2 p- 1 (x)dxds] + C p E[[ [ \Z^ x ' n - Z* s ' x \ 2 p- 1 {x)dxds] — » 0, 

Jt JR d Jt Js. d 

as n — ► oo. This will be used in the following theorem. 

Theorem 4.3 Under Conditions (H.1)-(H.4) ; if we define u(t,x) — Yl' x , where (Y*> x , Z\> x% ) is the 
solution of Eg. iTO]) , thenu(t,x) is the unique weak solution of Ea. ^.lb}) withu(T,x) = h{x). Moreover, 
u(s,X*' x ) = Y^ x , (cr*Vtt)(s,X*> x ) = Zl' x for a.e. s G [t,T], x G W l a.s.. 

Proof. From Proposition ^. 21 we only need to verify that this u is the unique weak solution of Eq. (|2.1Gp 
with u(T, x) = h(x). By Lemma [2~6l it is easy to see that (er*Vu)(t, x) = Z\ x for a.e. t G [0,T],i€ 
a.s.. Furthermore, by the generalized equivalence of norm principle again we have 



E[f I {\u{s,x)\ 2 + \{a*Vu){s,x)\ 2 )p- 1 {x)dxds\ 
Jo Jia d 

< C P E[ f T f (\u(s,X^ x )\ 2 + \(<T*Vu)(s,X^ x )\ 2 )p- 1 (x)dxd S ] 

JO JR d 

= C P E[ [ T f (\Y S °- X \ 2 + \Z° s > x \ 2 )p-\x)dxds] < oo. 
Jo Jm d 

Now we verify that u(t,x) satisfies (|2.17|) with u(T, x) = h(x) by passing the limit in L 2 (Q) to (|4.3p . 
We only show the convergence of the last term. The last term includes infinite dimensional integral, 
but 
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j=i Jt J Rd 

oo -t „ 

/ / 9j(s,x,u(s,x), (a*Vu)(s,x))W(s,x)dxd 1i f3 j {s)\ 2 ] 

j—lJt JR d 
n -t r 

<2£[|V / / ( gj (s,x,u n (s,x),(a*Wu n )(s,x)) - g J (s,x,u(s,x),(cr*Wu)(s,x)))>P(s,x)dxd ii /3 J (s)\ 
j=i J Rd 
oo r T r 

+2E[\ V / / g j (s,x,u(s,x),(a*Vu)(s,x))^(s,x)dxd^ j (s)\ 2 ] 
J=n+1 Jt Jr« 

< C p E[y2(Cj + ay) / / (|u n (i,x)-u(i,2:)| 2 + |((T*VM")(s ! a;)-(CT*Vu)(s,a;)| 2 )p- 1 (a;)(ia;(is] 

oo . 

+C P £[ | V / / (g j (s,x,u(s,x),(a*Vu)(s,x)) - g j (s,x,0,0))^(s,x)dxd^ j {s)\ 2 ] 
J=n+1 Jt Jr« 

oo fT r 

+C P E[\ J2 / / 9j (s,x,0, Q)^(s,x)dxd^ j (s) \ 2 }. 



j=n+l ' 



It is obvious that the first term tends to zero as n — * oo. The last two terms can be treated using a 
similar method as (|3.17j) and (|3.18|l . 

Therefore u(t,x) satisfies (|2.17p . so is a weak solution of Eq. (|2.16|) with u(T,x) — h(x). The 
uniqueness can be proved following a similar argument of Theorem 3.1 in Bally and Matoussi 3J. o 

5 Infinite horizon BDSDEs 

We consider the following BDSDE with infinite dimensional noise on infinite horizon, 

/oo />oo 
e- Kr f{r,X t r > x ,Y?< x ,Z t r > x )dr + J K e - Kr Y r tx dr 

poo poo 

- e- Kr g(r,X*' x ,Y?> x ,Z*' x )d*B r - e~ Kr (Z l r < x ,dW r ). (5.1) 

J S J S 

Here / : [0, oo) x I d x R 1 x R d — ► R 1 , g : [0, oo) x R d x M 1 x R d — > /^(K 1 ). Eq.(OJ> is equivalent to 

/oo poo 
e- Kr f(r,X t r ' x ,Y r t > x ,Z t r ' x )dr + J K C - Kr Y^ x dr 

°° pOQ pOQ 

-E/ B- Kr gj (r,X^X' x ,K' x )^^(r)- e- Kr (Z^ x ,dW r ). 
j=i J s J s 

We assume 

(H.5). Change ^ [0)T] to and t € [0, T] to i > in (H.2); 
(H.6). Change J Q T to J °° e - Ks in (H.3); 

(H.7). There exists a constant fx > with 2[i - K - 2C - Y^jLi c j > s -t- for an Y * > °> Yi,Y 2 € 
L 2 (M d ;M 1 ), I,Ze L 2 p (R d ;R d ), 
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(Y^x) -Y 2 (x))(f(t,X(x),Y 1 (x),Z(x)) - f(t,X(x),Y 2 (x),Z(x)))p-\x)dx 

I 

< - f i [ - Y 2 {x)\ 2 p~\x)dx. 



The main objective of this section is to prove 
Theorem 5.1 Under Conditions (H.4)-(H.7), Eg. \5.1\) has a unique solution. 
Proof. Uniqueness. Let (Y*> x , Z l s ' x ) and (Y*' x ,Zl' x ) be two solutions of Eq. (|5.1[) . Define 

y~'t.X \^t.X \rt,X iyt.X ryt.X 'yt^X 

s — J s — 1 s ' s — A s s ' 

f(s, x) = f( s , xt>*, £<•*) - /(«, x***, r^, 

ar) = <?( S , XI'", Y*'*, Z%*) - </;.-. X<>*, Y^ x , Z%*\ s > 0. 
Then for s > and a.e. x € M d , (Y^ x ,Zl' x ) and (Y/'*, £*> x ) satisfy 

/ dS?-« = -/>, x)ds + £~ fc- («, ar)^4-(«) + (^' x , <*W S ) 



liniT >oo e~ ^Yp' 111 = a.s. 



For a.e. x £ M. d , applying Ito's formula for infinite dimensional noise to e • K " s |y^' a! | , and by Young 
inequality and Conditions (H.5), (H.7), we obtain 

E[ e- Ks \Y*' x \ 2 p- 1 {x)dx]+{2fi-K-2C-y]C j )E[ / C - Kr \Y r ^ x f p- 1 (x)dxdr] 
JR d . =1 Js Jm d 

+ d-f2^)E[[ T I e- Kr \Z t - x \ 2 p-\x)dxdr] 



"2 



< E[ eT K1 \Y%*\ p- l {x)dx]. (5.2) 

Taking K 1 > K s.t. 2/i — if' — 2C — > as well, we can see that (|5.2[) remains true with K 

replaced by K' . In particular, 



E[ e- K ' s \Y s t ' x fp- 1 (x)dx}<E[ e- K ' T \Y^ x fp- 1 (x)dx}. 

JWL d JS. d 

Therefore, we have 

E[f e- K ' s \Y t s ' x \ 2 p- 1 {x)dx\<e^ K '- K ^ T E[[ e~ KT \Yp x \ 2 p- 1 {x)dx}. (5.3) 
Jm d JR d 

Since Y/' x , Y*> x e S 2 '~ K f] M 2 '- K ([0, 00); L 2 (R d ; M 1 )), so 

sup^[ / e- KT \Y*' x fp- 1 (x)dx] < £[sup [ e~ KT (2\Y*' x f + 2\Y*' x f)p- 1 {x)dx] < 00. 
Therefore, taking the limit as T — > 00 in (|5.3[) . we have 

E[f e- K ' s \Y^\ 2 p-\x)dx}=0. 

JR d 

Then the uniqueness is proved. 
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Existence . For each n G N, we define a sequence of BDSDEs (|3.1|) with h = and T = n and denote 
it by Eg. (13.11 ). It is easy to verify that for each n, these BDSDEs satisfy conditions of Theorem 13.21 



Therefore, for each n, there exists a (Y s *> x ' n , Z*> x >") G S 2 '°([0, n];L 2 p (R d ; R 1 )) <g> M 2 <°([0, n];L 2 p (R d ; R d )) 
which is equivalent to the space S 2, ~ K ([0, n); 

L 2 (M d ;M 1 ))®M 2 - iC ([0,n];I,2(R d ;K d )) and (Y s *> x '«, Z*> x ' n ) is the unique solution of Eq.(Ok)- That 
is to say, for an arbitrary tp £ C%(M d ;R 1 ), (Y s t ' x ' n , Z l s x - n ) satisfies 

e- Ks Y*' x > n ip(x)dx = f I e- Kr f(r,X^ x ,Y^ x ^,Z^ x ^)ip(x)dxdr 

Js JM d 



Ke- Kr Yt' x ' n <p(x)dxdr - V / / 



e -Kr 



e- Kr Z t r ' x ' n ip(x)dx,dW r ) P-a.s.. (5.4) 

Let (Y/\Z t ") t> „ = (0,0), then (Y*> x > n t Z*>*> n ) € S 2: ~ K f] M 2 ~ K ([0, oo);L 2 (R d ; R 1 )) M 2 ~ K 
([0,oo);L 2 (R d ;R d )). We will prove (Y^*'™, Z*' x > n ) is a Cauchy sequence. For this, let (Y s *> x ' m , Z t s ' x ' m ) 
and (Y s i,x > n , Z t s ,x ' n ) be the solutions of Eq. (|3.lD and Eg. 1)3.11 ,) respectively. Without losing any gen- 
erality, assume that m > n, and define 

■\>-t,x,m,n \rt,x,m \rt,x,n ryt^x^ra^n r^t.x^m iyt,x,n 

1 S — 1 S — 1 S 1 — A S ~ ^S ! 

Im.ni „ \ t^^t^^ ryt.x.m\ £( n vt,x \st,x.n r?t,x,n\ 

§f n («, z) = < (i is. .VI". Y s << x ' m , ^ x >™) - .VI'. Y s *' x '", Z<> x '«), a > 0. 
Consider two cases: 

(i) When n < s < m, Y s *> x ' ro > n = Y s *' x > m . Since K' I ' m ,^ I '" i ) is the solution of Eq.jSlL), we have 
for any m G N, 

f dY/' x - m = -f(s, Xl> x , Y s 4 ' x ' ro , Zl> x > m )ds + J^JLi 9]( s > x t' x > Y s *> x ' m , Z\- X ^ m )d^ ^(s) + <2W S ) 
\ Y^ x ' m = for s G [0,m], a.e. x G R d , a.s.. 

Noting that E[f™ \\g(r, Xp x , Y r *' x ' m , Z*< x ' m )|j 2 ^ (R1) dr] < oo for a.e. x G M d , we can apply Ito's formula 
to e - Kr \Y^' x ' m \ 2 for a.e. x G R d , then taking integration in R d x we have 

e- Ks \Y s t ' x ' m \ 2 p~ 1 {x)dx 

oo oo „ m „ 

+ (2/i-iV-2C-y - (1 + VQ)e) / / e-^l^^lV 1 ^)^ 

Js Js* 

1 °° °° r rn r 

+ (o-E^-E^ £ )/ / e- Kr \Z t r ' x ' m \ 2 p-\x)dxdr 

<C P / e- Kr \f(r,X t r ' x ,0,0)\ 2 p- 1 (x)dxdr 
Js Jm d 

n °° 
e" Kr E I* ( r > < °> «) | V 1 (*)^ 

V / / 2e- Kr Y r t > x > m g J {r,X t r > x ,Y r t < x < m ,Z t r < x ' m )p- 1 (x)dxd^ j (r) 

2e^ Kr Y r t ' a; ' m Z^^ m ^ 1 (a;)dx,dW r ). (5.5) 
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Note that the constant e can be chosen to be sufficiently small s.t. all the terms on the left hand side 
of (|5.5p are positive. By (|5.5[) . as n, m — > oo we have 

e- Kr \Y^ x ' m \ 2 p- 1 (x)dxdr]+E[ / e- Kr \Z t r ' x ' m \ 2 p- 1 (x)dxdr] 

J d Jn Jm d 

noo 
e~ Kr (\f(r, Xp*, 0, 0)| 2 + V \ 9j (r, Xp x , 0, 0)\ 2 )p- 1 (x)dxdr] — > 0. (5.6) 

3=1 

Note that the right hand side of (|5.6p converges to follows from the generalized equivalence of norm 
principle. Also using the B-D-G inequality to deal with (|5.5p in the interval [n,m], by (15. 6|) . as n, 
m — ► oo we have 

E[ sup f e- Ks \Y s t ' x ' m \ 2 p- 1 dx] 

n<s<m JM d 

e~ Kr (\f(r, Xp x , 0, 0)| 2 + V \ 9j (r, Xp x , 0, Q)\ 2 )p-\x)dxdr] 

J ' d 3 = 1 

+C p E[f f e- Kr (\Y r t ' x ' m \ 2 + \Z t r ' x ' m \ 2 )p- 1 (x)dxdr]^0. (5.7) 

Jn JtL d 

(ii) When < s < n, 

pn 00 p n pn 

yt, x , m ,n = yt, X ,rn + f*>«(r,x)dr l 0?>, X)£ fa (r) - / (^ TO '",dW r ). 
./ s ■ i J s J s 



3 = 1 ■ 

Apply Ito's formula to e-* r |Y J ?> a, ' m > n | 2 for a.e. x € R d , then 

OC 



/ e- Ks \Y s t > x > m < n \ 2 p- 1 (x)dx + {2fi-K-2C- Vc,) / / e-^S^'^lV^aOdxdr 

i ( \ * „ "\ / / „ — Kr I —t.x.m, n\2 „— 1 / 



'2 

3 = 1 



p (x)dxdr 



</ e- Kn \Y^ m \ 2 p-\x)dx~Y, / 2e"^y r ^ m ^7'"(r,a ; )p- 1 ( a; )dxdt^( r ) 

2e- Kr Yt' x < m - n Z t r - x ' m ' n p- 1 (x)dx, dW r ). (5.8) 
Taking expectation on both sides of (|5.8p , as n, m — ► oo, using (|5.7[) , we have 



JM d Js JM d 

< C P E[ sup / e-^lY,*^'" 1 ! 2 ^ 1 ^)^] — ►O. 

n<s<m JR d 

Also by the B-D-G inequality, (|5.7[) . (|5.8[) and (|5.9[) , as n, m — > oo, we have 

£[ sup / e-^lF^^lV 1 ^)^] < Cp^I sup / e~*'|l?' s ' m | V 1 —>0. 

Q<s<nJR d n<s<mJR d 

Therefore taking a combination of cases (i) and (ii), as n, m — ► oo, we have 



(5.9) 
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£[sup / e- Ks \Y 6 t ' x ' m - n \ 2 p" 1 (x)dx] +E[f f e~ Kr \Y r t ^' m - n \ 2 p- 1 (x)dxdr} 

s>0 JR d JO JW d 

+E[f f e- Kr \Z^ x ' m ' n \ 2 p- 1 (x)dxdr] — ► 0. 

JO JR d 

That is to say (Yf' x > n ,Z* a > x ' n ) is a Cauchy sequence. Take (Y*' x ,Zf x ) as the limit of (Yf> x > n ,Z t s > x > n ) 
in the space S 2 ~ K f] M 2 '~ K {[0, oo);L 2 (R d ; K 1 )) M 2 ~ K ([0, oo);L 2 (R d ; R d )) and we will show that 
(Y^ x ,Zl' x ) is the solution of Ea.([5T|l. We only need to verify that for arbitrary tp S C^K^jK 1 ), 
(Y* ,x , Z\' x ) satisfies (|2.14D . where (12. f 41 ) means a more general form of (|2 . 14[) with / and gj also 
depending on r € [0, oo). Since (Y*> x ' n ,Zl' x,n ) satisfies Eq. (|5.4|) . so we verify that Eq. (|5.4[) converges 
to Eq. (|2.14l in L 2 (Q) term by term as n — ► oo. We only show the infinite dimensional stochastic 
integral term: 

^IlE/ / e- Kr g j (r,Xp x X' x ' n ,K' x ' n M^)dx^0 j (r) 

~[ Js JR d 

°° p CO p 

"£/ / e- Kr g j (r,Xp x X' x ,K' X Mx)dxd^ j (r)\ 2 ) 
. : j s JR d 



3=1 ' 



poo p 
• x Jn JR d 



)(p(x)dxd^j(r)\ 21 



We see that each term in the above formula tends to zero as n — ► oo since 



£UE / / e^feia^F^^^-s^r,^^^^^))^,)^ 

j =l Js JR d 

< C P E[ f f e- Kr (\Y^ x ' n - Yy x \ 2 + \Z* r ' x ' n - Z^ x \ 2 )p^ 1 {x)dxdr] — ► 0, as n -> oo, 
Jo Jm d 



and 



^[lE/ / e-^lr,!^,^,^)^.)^^^)] 2 ] 

J - =1 Jn JR d 

< C P E[ ^ f e- Kr (\Y r ^\ 2 + \Z^\ 2 )p- 1 (x)dxdr] 
Jn Jm d 

poo p °° 

+C p / / E e ~ Kr igj( r » g > 0) 0)\ 2 p~ 1 (x)dxdr — ► 0, as oo. 

Jn JR d . =1 

That is to say (T s t,:E , Z' ,;l: ) s>0 satisfies Eq. (|2.14l ). The proof of Theorem 15. II is completed. o 

By similar method as in the proof of existence part case (i) in Theorem 15. 1[ we have the following 
estimation: 

Proposition 5.2 Let (Y£> x > n i Z%' x ' n ) be the solution of Eo. l\'6.ll ,), then under the conditions of Theo- 
rem 1 5. 1\ 

sup E [sup f e- Ks \Y s t ' x ^(x)\ 2 p~ 1 (x)dx}+ sup E[f f e~ Kr \Y^ n (x)\ 2 p- 1 (x)dxdr] 

n s>0 Jl J n Jo JR d 

+ supE[[ f e- Kr \Z t r ' x ' n {x)\ 2 p- 1 {x)dxdr] < oo. 
n Jo JH d 
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6 The continuity of the solution of the infinite horizon BDSDEs as the 
solution of the corresponding SPDEs 



Now we study BDSDE (J57T3J), a simpler form of Eq.(f57 

Proof of Theorem \2.1(A Since conditions here are stronger than those in Theorem 15.11 so there 
exists a unique solution (Y*' x , Zl' x ). We only need to prove -E[sup s>0 L d e~ pKs \Y*> x \ p p~ 1 (x)dx] < 
oo. Let ipN.p(x) = x% I{o<x<N} + f-/V £ 2 _ (x — N)I{ x >n\. We apply generalized Ito's formula to 
e~P Kr <p Np (ip M (Yj:' x )) for a.e. x € K d to have the following estimation 



~ ^e-^^^M^ 4 ^))!^^-)! 2 ^*-! 2 ^ 

- pKr <p'N, P (Y^)) I { -M<Y^<M} \Zr X ?dr 

< e-P KT ip N , p (i) M (Y*< x )) + £ e-P Kr ip NtP ^ M (Y r ^))^Yrnf(^ X , Yf'*, Z^ x )dr 

/T oo 
e- pKr ^ N ^M{Y^ x ))l { _ M<Y ^ <M} J2 \9j(K' x > Zp*)\ 2 dr 

3 = 1 

•>» j=l 

oo -t 

' pKr ^ N ^ M {Yy x ))i, M {Yy x )zr,dw r ). (6.i) 



i-T 



Note that limy^oo e pKT <Pn.p(' 1 Pm(Yt X )) — 0, so after taking limit as T ~ * oo, we take the integration 
on 12 x R d . As (Y?'°, Z.*'') <E S 2 ~ K f] M 2 '~ K ([0, oo); L 2 (R d ; R 1 )) M 2 '~ K ([0, oo); L 2 (K rf ; R d )) and 
i'iv p^mK^II'/'mK' 1 ) i s bounded, we can use the stochastic Fubini theorem and all the stochastic 
integrals have zero expectation. Using Conditions (A.1)'-(A.4)', and taking the limit as M — > oo first, 
then the limit as N — > oo, by the monotone convergence theorem, we have 

1 j=l Z j=l J s J ^ d 

OO OO n oo /» 

+ ?-(2p - 3 - (2p - 2) V a, - (2p - 2) £ "^)^[ / / e^l^f "V^lV 1 (x)dxdr] 
<C P / |/(a:,0,0)|^ 1 (a:)^ + C p / V) 0, 0)\* p- 1 (x)dx] < oo. (6.2) 

jR d JR d - =1 

Note that the constant e can be chosen to be sufficiently small s.t. all the terms on the left hand side 
of (|6.2p are positive. Also by the B-D-G inequality, Cauchy- Schwartz inequality and Young inequality, 
from (|6.1[) we have 

E[sup f e- pKs \Y^ x fp- 1 (x)dx] 

s>0 JK d 
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p p oo 

<C p \f(x 1 0,0)\Pp- 1 (x)dx + C p y2\g J (x 1 0,0)\Pp- 1 (x)dx 

pOO p pOO p 

+C P E[ / e- pKr \Y*' x f~ 2 \Z t r ' x \ 2 p- 1 {x)dxdr]+C p E[ / e' pKr \Y^ x \ p p~\x)dxdr}. 

JO JR d JO JR d 

So by J6l]), Theorem [230] is proved. o 

We need to prove two lemmas before giving a proof of Theorem 12.111 
Lemma 6.1 Under Condition (A. 3)', for arbitrary T > 0, t, t' £ [0, T], 

E[f f e~ Kr \X^' x -X t r ' x \ p p- 1 {x)dxdr] <C p \t'-t\% a.s.. 
Jo Jm d 

Proof. It is not difficult to deduce from Lemma 4.5.6 in [17j . so we omit the proof. o 

Lemma 6.2 Under Conditions (A.1)'-(A.4)', for arbitrary T > 0, t,t' £ [0,T], let (Y/' x ) s > , 
(Y*' x ) s>0 be the solutions of Eg. iTO]) , then 

E[sup f e- pKs \Y?' x - Y s t ' x \ p p- 1 (x)dx] < C p \t' - 
s>o Jm d 



Proof. Let 



y' vt ,a: \rt.x rr ryt .x 7t^x 

s — *s ~ J s ) — A B — & s , 

g 3 is)^9 3 {X^ x ,Y^ x ,Z^ x )-g 3 (Xr^ X ,Zr), s>0. 



Then 



f dF s = -/(s)rfs + E~ i + (Ja, dW s ) 

\ lim T ^oo e~ KT Y T = for a.c. x £ R d a.s.. 

First note that from Theorem 12.101 we know E[sup s>0 f Rd e~ pKs \Y s \ p p~ 1 (x)dx] < oo. Applying Ito's 
formula to e~ pKr \Y r \ p for a.e. x £ M. d (we leave out procedure of localization as in (|6.ip for simplicity) 
and taking integration on M. d , we have 



- pKs \Y s \ p p- 1 (x)dx 

P(P~ 1) 



(pp — pK — pC — 



2 

j 



Vc r 3 £ ) / e- pKr \Y r \ p p- 1 (x)dxdr 

= 1 Js Jwi d 



< 



°Q pOO p 

+?(2p-3-(2p-2)V aj -) / / e-P^^M^lVV^dr 

4 5=1 J* J Rd 

C p I e- pKr \X r \ p p- 1 {x)dxdr -p^ / / e~ p,rr |r r |'- a ? r fl i (r)p- 1 (a!)d!i:<?4-W 

=1 Js Jm d 

-P f (f e- pKr \Y r \P- 2 Y r Z r p- 1 (x)dx,dW r ). (6.3) 



Note that the constant e can be chosen to be sufficiently small s.t. all the terms on the left hand side 
of (|6.3j) are positive. Taking integration on f2 on both sides of (|6.3|) . by Lemma l6~Tl we have 
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pOQ p pOO p 

E[ / e- pKr \Y r \ p p- 1 (x)dxdr] + E[ / e~ pKr \Y r \ p - 2 \Z r \ 2 p- 1 {x)dxdr] 

Js JM. d Js JR d 

POO P 

<C P E[ / e- pKr \X r \ p p- 1 (x)dxdr] 

Js JR d 

<C p \t'-t\%. (6.4) 
Also by the B-D-G inequality, from (|6.3p and (|6 .4|) . we have 

£[sup / e- pKs \Y s \ p p~ 1 (x)dx] 

s>0 jR d 

poo p poo p 

<C P E[ / e- pKr \X r \ p p-\x)dxdr]+C p E[ / e~ pKr \Y r \ p p- 1 (x)dxdr] 
Jo Jm d Jo Js. d 

+C P E[ f e- pKr \Y r \ p ' 2 \Z r \ 2 p- 1 (x)dxdr] 

Jo JR d 

< C p \t'-t\%. 

o 

Proof of Theorem \2.11\ By Lemma 16. 2\ we have 

E([snp f e ~ 2Ks \Y^ x - Y^ x \ 2 p- 1 (x)dx])i 

s>0 JR d 

<C p E[sup [ e- pKr \Yj'' x -Y*> x \ p p- 1 {x)dx]{ [ p- 1 (x)dx) E ^ 



s>0 

< C p \t'-t\%. 

Noting p > 2, by the Kolmogorov continuity theorem (see [T7]), we have t — > y*> x is a.s. continuous 
for t £ [0,T] under the norm (sup s>0 J Rd e~ 2Ks \ ■ \ 2 p~ 1 (x)dx) ' . Without losing any generality, assume 
that t' > t. Then we can see 

lim( / e- 2Kt '\Y^ x - Y t t / x \ 2 p- 1 {x)dx)^ < lim(sup / e - 2Ks \Yf- x - Y*' x \ 2 p- 1 (x)dx)^ = a.s.. 

* ,_>t JR d *'^ t s>0 jR d 

Notice t' € [0,T], so 

lim(/ ly^^-^lV 1 ^)^)* =° a - s - (6.5) 

Since Y*'' G S 2 "^([0, oo); K 1 )), if,'' is continuous w.r.t. f in L^K^K 1 ). T h a t is to say for 

each t, 



lim( / \Y*<* - Y t t ' x \ 2 p~ 1 (x)dx)i = a.s.. (6.6) 



Now by d63]) and [[61) 



iim(/ ir/^-r^iV 1 ^)^)^ 

4 JR d 

< hm ( / |y/'» - Y*>*\*p- X (x)dx)$ + lim( / \Y$* - Y^ x \ 2 p-\x)dxY 
= a.s.. 



Stationary Solution of SPDEs 



35 



For arbitrary T > 0, < t < T, define u(t, •) = F t *'', then u(t, •) is a.s. continuous w.r.t. t in L 2 p (R d ; M 1 ). 

Since K*'' G <S 2 '~ K ([0, oo); Lp(R d ; M 1 )), is^U®^" measurable and E[J Rd \Y% ' x \ 2 p~ l [x)dx] < 
oo. It follows that Condition (H.l) is satisfied. Moreover, Conditions (A.1)'-(A.3)' are stronger than 
Conditions (H.2)-(H.4), so by Theorem |43I u(t, x) is a weak solution of Eq. Q2.16p . Theorem f2. Ill is 
proved. o 
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